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ABSTRACT 

We study zCOSMOS-bright optical spectra for 609 Spitzer/MIPS 24 /im -selected galaxies with 
SW^m > 0.30 mJy and / < 22.5 (AB mag) over 1.5 deg 2 of the COSMOS field. From emission- 
line diagnostics we find that: 1) star-formation rates (SFR) derived from the observed Ha A6563and 
Hf3 A4861 lines underestimate, on average, the total SFR by factors ~ 5 and 10, respectively; 2) 
both the Calzetti et al. and the Milky Way reddening laws are suitable to describe the extinction 
observed in infrared (IR) sources in most cases; 3) some IR galaxies at z < 0.3 have low abundances, 
but many others with similar IR luminosities and redshifts are chemically enriched; 4) The average 
[OIII] A5007 /H/3 A4861 ratios of vLf^ > 10 11 L© galaxies at 0.6 < z < 0.7 are - 0.6 dex higher than 
the average ratio of all zCOSMOS galaxies at similar redshifts. Massive star formation and active 
galactic nuclei (AGN) could simultaneously be present in those galaxies with the highest ionising 
fluxes; 5) ~ 1/3 of the galaxies with metallicity measurements at 0.5 < z < 0.7 lie below the general 
mass-metallicity relation at the corresponding redshifts. The strengths of the 4000 A break and the 
Hi5 EW of our galaxies show that secondary bursts of star formation are needed to explain the spectral 
properties of most IR sources. The LIRG and ULIRG phases occur, on average, between 10 7 and 10 s 
years after the onset of a starburst on top of underlying older stellar populations. These results are 
valid for galaxies of different IR luminosities at 0.6 < z < 1.0 and seem independent of the mechanisms 
triggering star formation. 

Subject headings: infrared: galaxies - galaxies: abundances galaxies: active - galaxies: evolution 
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Spectroscopic surveys are probably the finest diagnos- 
tics to probe the nature of galaxies. They provide ac- 
curate redshift determinations and directly reveal the 
physical and chemical processes governing galaxy evolu- 
tion through cosmic time. The high cost of spectroscopic 
campaigns has prevented for many years the follow-up of 
a significant number of sources in large areas of the sky. 
Yet, extensive spectroscopic surveys are necessary in or- 
der to characterise in detail different galaxy populations 
in a representative way. 

Among galaxy populations, infrared (IR) galaxies have 
been a subject of major interest in observational cosmol- 
ogy since long before the discovery of the extragalactic 
IR background (Puget et al. 1996). The most recent de- 
terminations indicate that the integrated emission from 
IR galaxies has a significant contribution to the extra- 
galactic background light, comparable to that produced 
at optical wavelengths (Dole et al. 2006). As signposts of 
star-formation and active galactic nucleus (AGN) activ- 
ity, IR galaxies play a major role in reconstructing galaxy 
formation, evolution and the history of stellar mass as- 
sembly. 

Surveys conducted with current-generation IR facil- 
ities -such as the Spitzer Space Telescope (Werner et 
al. 2004) and now also the Akari Telescope (Matsuhara 
et al. 2006)- are rapidly improving our understanding of 
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the nature and evolution of IR galaxies up to high red- 
shifts. It is now well established that luminous and ultra- 
luminous IR galaxies (LIRGs and ULIRGs; Sanders & 
Mirabel 1996) are much more common at high (z > 0.5) 
rcdshifts than in the local Universe and the star forma- 
tion rate (SFR) density strongly increases up to red- 
shift z ~ 1 (e.g. Le Floc'h et al. 2005; Caputi et 
al. 2006a, 2007). These results have been previously 
suggested by studies conducted with the Infrared Space 
Observatory (ISO) (e.g. Aussel et al. 1999; Chary & 
Elbaz 2001; Franceschini et al. 2001). In addition, multi- 
wavelength photometric observations of IR galaxies have 
produced further information on their different proper- 
ties: already-assembled stellar masses (Daddi et al. 2005; 
Caputi et al. 2006a, b; Papovich et al. 2006), character- 
istic spectral energy distributions (SEDs; e.g. Dale et 
al. 2005; Rowan-Robinson et al. 2005; Rocca-Volmerange 
et al. 2007; Takagi et al. 2007; Bavouzet et al. 2008), and 
the AGN fraction among IR sources (e.g. Alonso-Herrero 
et al. 2006; Treister et al. 2006). 

Spectroscopic data are necessary to address some other 
questions, such as e.g. the chemical composition and 
metallicities of IR galaxies at different redshifts. Also, 
they are useful to understand the distribution of dust 
within a galaxy, by comparing the degree of reddening 
independently inferred from the IR flux and from the 
spectral line decrement. Finally, spectroscopic data can 
help to disentangle the presence of AGN among IR galax- 
ies. 

Some of the first spectroscopic surveys of Infrared 
Astronomical Satellite (IRAS) sources studied the nu- 
clear regions of nearby LIRGs and ULIRGs (Kim et 
al. 1995; Veilleux et al. 1995, 1999). Later, Franceschini 
et al. (2003); Flores et al. (2004) and Liang et al. (2004) 
studied the spectra of /S'O-selected IR galaxies at z ~ 
0.6 — 0.7. More recently, Choi et al. (2006) analysed 
spectroscopic star formation rates (SFR) and extinction 
properties of Spitzer sources with 24 fim detections. Pa- 
povich et al. (2006b) obtained spectroscopic redshifts for 
a very large sample of 24 /zm -selected galaxies, but they 
presented a very limited analysis of the corresponding 
spectra. Even though all these works investigated dif- 
ferent spectral properties of IR sources of different lumi- 
nosities and redshifts, none of them provides a complete 
analysis of a large homogeneous sample of IR-selected 
galaxies from redshifts z ~ to z ~ 1. This is the goal 
of this paper. 

The Cosmic Evolution Survey (COSMOS; Scoville et 
al. 2007) has been designed to probe galaxy evolution, 
star formation and the effects of large-scale structure up 
to high redshifts. The COSMOS field is defined by its 
Hubble Space Telescope/ Advanced Camera for Surveys 
(EST / AGS) 2 deg 2 coverage. There are multiple follow 
up observations carried out in the COSMOS field, rang- 
ing from X-rays to radio wavelengths. Among the photo- 
metric imaging observations, COSMOS includes the full 
and homogeneous coverage of the field with the Spitzer 
Infrared Array Camera (Spitzer /IR AC; Fazio et al. 2004) 
and the Multiband Photometer for Spitzer (MIPS; Rieke 
et al. 2004), as part of Spitzer cycle-2 and 3 Legacy Pro- 
grams (Sanders et al. 2007). 

COSMOS comprises also a large spectroscopic follow- 
up program being performed with the Visible Multiob- 
ject Spectrograph (VIMOS) on the Very Large Telescope 



(VLT): zCOSMOS (Lilly et al. 2007). The zCOSMOS 
survey consists of two parts: 1) the observation of an 
/ < 22.5 AB mag limited sample of 20,000 galaxies over 
1.7 deg 2 of the COSMOS field (zCOSMOS-bright) and 2) 
a sample of 10,000 galaxies in the central 1 deg 2 , colour- 
selected to have redshifts 1.4 < z < 3.0 (zCOSMOS- 
deep). 

This paper constitutes the first of a series of studies 
we are carrying out to analyse the optical spectral prop- 
erties of 24 /jm -selected galaxies in the COSMOS field 
up to high redshifts. Here, we analyse the sample of 611 
24 fim sources with S , 24 Atm > 0.30 mJy that have been 
observed so far with zCOSMOS-bright among a total of 
~ 10, 000 (10k) sources. This is one of the largest spec- 
troscopic samples of mid-IR-selected galaxies analysed to 
date in the redshift ranges < z < 1 and < z < 3 
for normal galaxies and AGN, respectively. The lay- 
out of this paper is as follows: in Section Sj2j we de- 
scribe in detail our datasets and the cross-correlation of 
24 /im and zCOSMOS-bright sources. In Section £j3l we 
compute rest-frame mid-IR luminosities for our galaxies. 
We present our results on the analysis of the spectra in 
the following sections: mean and typical dispersion of 
the spectra in different redshift and IR luminosity bins 
O, comparison of different SFR indicators (§5j), the 
study of dust properties of our sample (<j6]) and line ra- 
tio diagnostics and derived metallicities (jFJ. Later on, 
in Section ^Hl we analyse the AGN present in our sam- 
ple. In section $9l we show how the optical spectra can 
put constraints on the average star formation histories of 
IR galaxies and, in section M 101 we analyse galaxy can- 
didates for being at the earliest stages of a burst of star 
formation. Finally, we summarise and discuss our results 
in Section fTTI We adopt throughout a cosmology with 
H = 70kms- 1 Mpc-\ Q M = 0.3 and Q A = 0.7. 

2. DATASETS 

2.1. The zCOSMOS-bright 10k sample 

The zCOSMOS survey (Lilly et al. 2007) is a very large 
spectroscopic program being performed with VIMOS (Lc 
Fevre et al. 2003) on the VLT. VIMOS is a multi-slit 
spectrograph with four non-contiguous quadrants which 
cover 7x8 arcmin 2 each. zCOSMOS is designed to have 
a uniform pattern of pointings and a multiple-pass strat- 
egy, which guarantee a uniform coverage of the field. We 
refer the reader to Lilly et al. (2007) for further technical 
details about this survey. 

In zCOSMOS-bright, the targets are selected randomly 
from a complete I < 22.5 AB mag catalogue (the 'parent 
catalogue'). This part of the program includes also a set 
of fainter targets which have been selected for specific 
reasons, e.g. X-ray or radio sources. All the observations 
are performed with the R ~ 600 VIMOS MR grism, with 
a spectral coverage of (5500 — 9500) A. 

At the moment of writing, roughly one half of the 
zCOSMOS-bright targets have already been observed 
and their spectra have been reduced using the VIPGI 
software (Scodeggio et al. 2005). Although the data re- 
duction is mostly an automated process, the redshift de- 
termination of each source is checked manually on an 
individual basis by two independent reducers and the 
results are reconciled later on. The resulting redshift 
catalogue and set of reduced spectra for this first com- 
pleted part of zCOSMOS-bright contains 10,644 sources 



The optical spectra of 24 /jm galaxies in COSMOS: I 



3 



over 1.5 deg 2 ('the 10k sample' hereafter). 10,580 out 
of these sources have / < 22.5 AB mag (Lilly et al., in 
preparation) . 

The confidence in the redshift determination of each 
source is qualified with a flag, whose values can be: 4 
(completely secure redshift); 3 (very secure redshift, but 
with a very marginal possibility of error); 2 (a likely 
redshift, but with a significant possibility of error); 1 
(possible redshift); (no redshift determination); and 9 
(redshift based on a single secure narrow line, which is 
usually [Oil] A3727or HaA6563). The addition of +10 
to the flag indicates a broad- line AGN (BLAGN). 

In addition, photometric redshifts for the entire sam- 
ple have been obtained with The Zurich Extragalac- 
tic Bayesian Redshift Analyzer (ZEBRA; Feldmann et 
al. 2006), using the multiple UV through near-IR wave- 
band data available for the COSMOS field. Tests per- 
formed on objects with duplicate spectra of different 
quality show that good-quality photometric redshifts as 
those obtained with the COSMOS datasets can be use- 
ful to confirm less-secure spectroscopic redshifts, such as 
those with flag=2, 1 or 9 (see Lilly et al. 2007). 

2.2. The SCOSMOS MIPS shallow 24 /im catalogue 

The SCOSMOS survey (Sanders et al. 2007) com- 
prises IRAC 3.6, 4.5, 5.8 and 8.0 /im and MIPS 24, 70 
and 160 /im observations including the entire 2 deg 2 of 
the COSMOS-ACS field, as part of the Spitzer cycle- 
2 and 3 Legacy Programs. In cycle-2, the COSMOS 
field has been mapped at 24 /on down to a completeness 
flux density S^m = 0.30 mJy (SCOSMOS-shallow). A 
smaller 30x20 arcmin 2 region has been observed down 
to a flux density limit #24 ~ 0.08 mJy as a verifica- 
tion field. The observations performed in cycle-3 extend 
these deeper maps to the entire COSMOS area. 

The source extraction on the 24 /im maps has been per- 
formed with the IDL version of the DAOPHOT pack- 
age. A point spread function (PSF) fitting technique for 
measuring the photometry has been necessary to deal 
with blending problems. Further details on the SCOS- 
MOS survey as well as the data reduction and source 
detection procedures are given by Sanders et al. (2007). 
The final SCOSMOS-shallow 24 /im catalogue contains 
9,807 sources with S^i^m > 0.30 mJy over an area 
1.75x1.97 deg 2 . Out of them, 3,150 sources lie within 
the zCOSMOS-bright lOk-sample field (see Figure [J). 

We explain in detail the selection effect produced by 
the zCOSMOS-bright 10k sample in next section and 
summarise the numbers of sources in Table [TJ 

2.3. Cross- correlation of the catalogues 

We firstly checked the fraction of 24 /im sources present 
in the zCOSMOS-bright parent catalogue, i.e. the cat- 
alogue from which the spectroscopic targets are taken 
randomly. This is necessary to assess how representative 
the zCOSMOS-bright 10k sample is for the identification 
of sources in SCOSMOS-shallow. 

We cross-correlated the SCOSMOS-shallow 
24 /im catalogue with the zCOSMOS-bright parent 
catalogue, using a matching radius of 2 arcsec. We 
found that 2,084 out of 3,150 (-66%) 24 /mi sources 
have a counterpart in the zCOSMOS-bright parent 
catalogue limited to / < 22.5 AB mag. Four additional 



sources have fainter /-band counterparts which are also 
in the potential list of targets of zCOSMOS-bright. The 
24 /im sources with associations in the zCOSMOS-bright 
parent catalogue have a similar 24 /im flux distribution 
as the entire sample of 3,150 24 /im sources lying in the 
zCOSMOS-bright lOk-sample field (see Figure HJ. Thus, 
the 24 /im sources identified with zCOSMOS-bright 
objects can be considered as representative of the 
majority (~ 66%) of the S^/un > 0.30 mJy population 
at different fluxes. 

The 24 /im sources without an association in the 
zCOSMOS-bright parent catalogue have I > 22.5 AB 
mag. Photometric redshifts obtained with ZEBRA for 
I < 24 AB mag sources (Oesch et al., in preparation) in- 
dicate that the majority ( > 80%) of the MlPS-shallow 
sources not present in the zCOSMOS-bright parent cat- 
alogue are at z > 1. 

703 out of the 3,150 24 /im sources in the zCOSMOS- 
bright field have been observed so far and form part of 
the zCOSMOS-bright 10k sample. 

Multiple associations are a relatively minor problem 
in the identification of 24 /im sources. To verify the reli- 
ability of each zCOSMOS-bright counterpart, we cross- 
correlated the 24 /im catalogue with a general / < 25 AB 
mag catalogue available for the COSMOS field (Capak 
et al. 2007). In 35 out of 703 cases (< 5%), we found 
that there is an / < 25 mag source which is closer to 
the 24 /im centroid than the zCOSMOS-bright associated 
source. We assumed that the 24/im-zCOSMOS associa- 
tion was false in these cases, and we have not further con- 
sidered these sources in our analysis. In other 73 out of 
703 cases (~ 10%), there is an / < 25 mag source which 
is also within 2 arcsec distance of the 24 /im centroid, 
but is further away than the zCOSMOS counterpart. 
We considered these 24 /zm -zCOSMOS associations as 
'likely'. The inclusion of these sources does not in- 
troduce any significant change in our results. Finally, 
595 out of 703 (~ 85%) are secure one-to-one 24 /zm- 
zCOSMOS associations. The resulting sampling rate for 
the 24 //m galaxies with / < 22.5 AB mag is — 1/3, sim- 
ilar to the overall sampling rate of the parent catalogue 
achieved so far with the total 10k sample. 

Having 668 24 /im sources with secure or likely 
zCOSMOS-bright-lOk associations, we had to determine 
how many of them had good-quality zCOSMOS spectro- 
scopic redshifts and spectra. We found that 482 out of 
these 668 sources have a zCOSMOS redshift with qual- 
ity flags 3 or 4. Additional 129 sources have less secure 
spectroscopic redshifts (flags=2,l or 9), but which are 
confirmed by ZEBRA photometric redshifts. The analy- 
sis of a subset of repeated zCOSMOS spectra of different 
quality has showed that > 90% of the spectroscopic red- 
shifts with flag=2,l or 9 which are confirmed by ZEBRA 
are correct (Lilly et al., in preparation). 

In total, 611 out of the 668 (i.e. > 90%) considered 
associations have secure redshifts and reasonable-quality 
spectra to perform line measurements. Their 24 /im -flux 
distribution is similar to that of all galaxies identified 
in the zCOSMOS-bright parent catalogue (compare thin 
and thick- line shaded histograms in Figure [2]) . 2 out of 
the 611 sources with secure redshifts have been identified 
with galactic stars. The remaining 609 galaxies consti- 
tute the sample analysed in this work. 

The photometric redshifts obtained with ZEBRA for 



4 



K. I. Caputi et al. 



the 57 spectroscopic failures suggest that 70 and 30% of 
them are at z < 1 and z > 1, respectively. 

3. REST-FRAME 24 LUMINOSITIES 

We classify the 609 24 /im galaxies in our spectro- 
scopic sample according to their rest-frame IR lu- 
minosities and zCOSMOS redshifts. We computed 
rest- frame 24 /j,m luminosities vL 2 ^ Aim using the observed 
24 /im fluxes, the zCOSMOS redshifts and model-derived 
k-corrections; i.e. vL 24t * m = ^47rfc(A)524 ^111^(2:), where 
fc(A) is the k-correction, 524^m is the 24 /im flux and 
dh{z) is the luminosity distance. We show the rest- 
frame 24 /im luminosities versus zCOSMOS redshifts for 
our galaxy sample in Figure [3] 

We used the Lagache et al. (2003, 2004) IR SED mod- 
els to compute the mid-IR k-corrections. As most of our 
galaxies lie at redshifts z < 1, the use of other models 
(e.g. Chary & Elbaz 2001; Dale & Helou 2002) would 
produce rest-frame 24 /im luminosities which are consis- 
tent within a factor 2. The discrepancies can be consid- 
erably larger at higher redshifts (Caputi et al. 2007). 

Traditionally, it is preferred to classify IR galaxies ac- 
cording to their bolometric IR (5-1000 /im) luminos- 
ity -L{^j . However, we decided to classify our sam- 
ple using their rest-frame 24 /im luminosities vL 2A Mm , 
as the monochromatic 19 IR luminosity is a less model- 
dependent quantity. As we explain below, we do apply 
recipes to convert the rest-frame vL 2A Mm to bolometric 
IR luminosities when necessary, e.g. for the computation 
of star formation rates (SFR). Roughly, a factor ~ 10 
should be applied to convert vL 2i Mm into L^ L (although 
this factor decreases with increasing i/L 24flm ). Thus, 
galaxies with vL 2A ^ m < 10 10 L Q , 10 10 < vLf»™ < 
10 11 L Q and vL 24fim > 10 11 L Q would roughly corre- 
spond to IR normal galaxies, LIRGs and ULIRGs, re- 
spectively. 

The Lagache et al. and other SED models are applica- 
ble to IR galaxies dominated by star formation. The SED 
of some AGN, instead, have a power-law shape f v oc v a 
(a < 0; see e.g. Elvis et al. 1994; Alonso-Herrero et 
al. 2006). 

We investigated the presence of IR power-law SED 
galaxies within our sample, by analysing the correspond- 
ing Spitzer/IRAC photometry. We classified as IRAC 
power-law sources those galaxies having fluxes in the four 
IRA C channels 3.6, 4.5, 5.8 and 8.0 /zm that are consis- 
tent with a single power law /„ oc v a (a < 0), within the 
error bars. For this, we considered a minimum error bar 
of 0.10 mag in each IRAC channel. 

45 out of the 609 24 /im galaxies in our sample have 
IRAC power-law SED (see Figure EJ). 39 out of the 45 
IRA C power-law SED sources are detected in the XMM- 
Newton X-rays maps of the COSMOS field (Hasinger et 
al. 2007) and 30 are BLAGN, as classified from their 
zCOSMOS spectra. We further discuss the IRAC power- 
law sources and other AGN present in our sample in Sec- 
tion Sj8j 

4. THE SPECTRA 

Before performing different spectral line diagnos- 
tics, we did some general inspection of the qualita- 

19 In our case, the rest-frame IR luminosities are convolved with 
the 24 fim filter response function rather than monochromatic. 



tive and quantitative spectral characteristics of our 
24 /im galaxies. 

Each individual zCOSMOS spectrum covers the ob- 
served wavelength range (5500 - 9500) A and has a dis- 
persion of 2.55 A. The quality of the spectra is quite 
uniform up to wavelengths A ~ 8000A. At longer wave- 
lengths, the quality is degraded by fringing, which can 
become quite severe in some cases beyond A - 8500A. 
As we explain later, we only perform line measurements 
on individual spectra at A < 8500A, and all the mea- 
surements between A = 8000 and A = 8500 A have been 
manually checked on an individual basis. 

Stacked spectra can give a good qualitative idea of 
the mean optical spectra of 24 /im galaxies. Figure HJ 
[5] and [6] show the rest-frame average stacked spectra 
of our 24 /im galaxies with vL 24 ^ < 10 10 L Q , 10 10 < 
vL 2i^m < i O "L0and vLf^ m > 10 11 L Q , respectively, 
in different redshift bins. We constructed the stacked 
spectra on sets of 38 to 89 galaxies, depending on the 
redshift and IR luminosity bin. The zCOSMOS BLAGN 
have been excluded for the stacking. The wavelength res- 
olution in all the stacked spectra is 2A (rest-frame). We 
only show the spectra up to rest-frame wavelengths cor- 
responding to observed A < 8500 A, i.e. we excluded the 
regions most affected by fringing. However, as we show 
in Section [J71 even the longest wavelength regions show 
reasonably good average spectra when sufficient numbers 
of sources are stacked. This is due to the fact that the 
fringing pattern adds incoherently for different sources, 
so no systematic noise is propagated into the average 
stacked spectra. 

We obtained the average spectrum in each bin by re- 
normalising the individual rest-frame spectrum of each 
source to the average value of a featureless region of the 
continuum (which was chosen depending on the redshift 
bin). In this way, all the individual spectra in a given 
redshift and IR luminosity bin are put on a same scale 
before stacking. We smoothed out regions lying on top 
of the main atmospheric absorption lines, except when 
source emission lines were present. At each rest- frame 
wavelength, we excluded the 5% smallest and largest val- 
ues before computing the average. This sigma-clipping 
procedure helps to clean the stacks for possible remaining 
spurious lines in the individual spectra. 

These average spectra of 24 /im galaxies in different 
redshift and IR luminosity bins show that: 

• as expected, all emission lines characterising star- 
forming galaxies are present. We note that these 
emission lines are a property of the average spectra 
of IR galaxies, but we do not necessarily observe all 
these lines in every spectrum on an individual basis 
(because e.g. they are much extincted by reddening 
in some cases). On the average spectra, we also see 
some absorption lines characteristic of old stellar 
populations, as NaD or Call H & K. This means 
that different generations of stars are present in 
many IR galaxies. 

• the average line ratios vary as a function of IR lumi- 
nosity and redshift. We further analyse this point 
in Section EjTj 

• high-order Balmer absorption lines are clearly 
present in 24 /im galaxies. These lines are produced 
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by short-living ( < 1 Gyr) A-type stars. This in- 
dicates that 24 /im galaxies not only are instanta- 
neously forming stars, but have been also forming 
stars for some time during the last Gyr. We explore 
this issue in more detail in Section S}9] 

Although the stacked spectra give a good qualitative 
idea of the average spectral properties of 24 /im galaxies 
of different IR luminosities and redshifts, they do not 
contain information about the variety of strengths of the 
spectral features among 24 /im galaxies in a same redshift 
and luminosity bin. 

We measured the fluxes and equivalent widths (EW) 
of emission lines in our zCOSMOS spectra by direct in- 
tegration on the rest-frame spectra. We do not measure 
emission lines lying beyond observer-frame A = 8500 A, 
to avoid being severely affected by fringing. For lines at 
observer-frame 8000 < A < 8500 A, we manually checked 
the measurements on an individual basis using the IRAF 
package 'splot'. Wc also checked manually those lines 
lying on top of or near strong atmospheric absorption 
features, e.g. [OII]A3727at z ~ 0.7. In all the follow- 
ing, except when otherwise stated, we only present and 
analyse line measurements when the corresponding rest- 
frame EW are > 5 A, as those are the lines with the most 
secure measurements. 

The error bars have been computed as the sum of two 
components. The first error-component has been ob- 
tained by considering the variations of the continuum in 
a narrow region around the line and setting this to a min- 
imum of 5%. On the other hand, ~ 10% of the sources 
on which we performed line measurements have dupli- 
cate good-quality spectra. We used these repeated line 
measurements to work out a second component for the 
error bars. The error bars computed in the first step have 
been then systematically shifted, in such a way that their 
median values coincided with the medians derived from 
the analysis of duplicates. The systematic percentual er- 
rors we added to the line fluxes are 12, 15, 12, 9 and 12% 
for [Oil] A3727 , H/3 A4861 , [OIII] A5007 , Ha A6563 and 
[Nil] A6584 , respectively. For the EW measurements, 
the systematic percentual errors added are 17, 20, 17, 10 
and 18%. 

All Balmer emission lines have been corrected for stel- 
lar absorption. To do this, we fitted the continuum 
of each spectrum with synthetic stellar SED from the 
Bruzual & Chariot 2007 library (Bruzual & Chariot 2003; 
Bruzual 2007). The medians of the stellar absorption cor- 
rections for HaA6563and H/3A4861in our galaxies are, 
respectively, (1.8 ± 0.5) A and (3.1 ± 0.3) A. 

In addition, we applied aperture corrections to all the 
spectral lines, as the slits in the VIMOS masks have a 
width of 1 arcsec. We convolved each zCOSMOS spec- 
trum in our sample with the SUBARU R and /-band 
filters and compared the resulting magnitudes with the 
total magnitudes of our galaxies, as available in the COS- 
MOS optical photometry catalogues (Capak et al. 2007). 
The difference between the two (spectrum-derived and 
total) magnitudes gave us the aperture correction in each 
case. For the majority (> 75%) of our galaxies, the 
optical-flux correction factors are between 1 and 4. For 
less than 2% of the sample, in all cases z < 0.3 galaxies, 
the correction factors can be as large as a factor 10 to 
20. 



Figure [7| shows the EW of typical emission lines in 
the spectra of 24 fim galaxies as a function of redshift. 
Only cases with EW > 5 A are presented in this fig- 
ure. Sources with Ha A6563 EW > 5 A are more than 
80% of all our sources at z < 0.3 (the redshift bin 
in which we performed Ha A6563 measurements on an 
individual basis). Instead, sources with H/3A4861and 
[OIII] A5007EW > 5 A constitute ~ 35 to 40% of all our 
sources at 0.2 < z < 0.7. 

From figure \7\ we see that the spectral lines of 
24 /im galaxies display a large range of EW values. Even 
galaxies with comparable IR luminosities at similar red- 
shifts have lines whose EW differ, in some cases, by a 
factor 5 or more. 

There is a wide range of possible Ha A6563 EW values 
for galaxies with vl?^^ < 10 10 L Q up to z = 0.3. This 
is probably a spatial sampling effect, because the VIMOS 
slit only samples a small region of 1 arcsec width. Within 
a star-forming galaxy, different regions such as HII nu- 
clei, hotspots or disk HII regions are characterised by 
different Ha A6563 EW distributions (e.g. Kennicutt et 
al. 1989). At low redshifts, these different regions can be 
resolved and they can individually dominate the optical 
light collected in the 1 arcsec-wide slit. In contrast, given 
the size of the MIPS PSF at 24 /im (~ 6 arcsec), the IR 
luminosity corresponds to the IR light in the integrated 
galaxy. The large dispersion in EW is not seen anymore 
for H/3A4861at higher redshifts, where the different re- 
gions cannot be resolved (see Figure [7]). However, the 
[OIII] A5007 EW do span a wide range of values at these 
higher redshifts, indicating the existence and different 
importance of ionising sources in IR galaxies. 

5. COMPARISON OF STAR FORMATION RATE 
INDICATORS 

When all the ultra-violet (UV) light produced by young 
stars in a star-forming galaxy is absorbed and re-emitted 
by the surrounding dust, the IR bolometric (5 — 1000 /im) 
luminosity Lj^i is directly related to the SFR (Kennicutt 
1998). Sometimes, however, part of the light produced by 
young stars is not absorbed and can directly be observed 
at UV wavelengths (for a discussion see e.g., Buat et 
al. 2007). This mainly depends on the geometry and 
mixing of the dust with young stars. Thus, an optimal 
estimator of SFR should take into account both the IR 
and UV contributions. 

Our aim is to compare the SFR derived from emission 
lines with no dust-correction in the zCOSMOS spectra, 
with the fiducial total (IR+UV) SFR. This allows us, 
in the next section, to quantify the degree of extinction 
characterising 24 /im galaxies. 

We consider here that the total SFR of our 
24 /im galaxies is given by the sum of the IR and UV- 
derived SFR: 

S , Fi? m (M yr" 1 ) = 1.72 x RU 10 x , (1) 

SFRuv(M Q yr- 1 ) = lA x 10~ 28 x L^ v ; (2) 

where Lj^j is the bolometric IR luminosity expressed in 
L and is the UV luminosity expressed in erg s _1 
Hz -1 , which is approximately constant between 1500 and 
2800 A for a Salpeter (1955) initial mass function (IMF; 
see Kennicutt 1998). Both equations assume a Salpeter 
IMF over stellar masses (0.1-100) M . 
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We computed the bolometric IR luminosity of our 
galaxies based on the rest-frame 24 /im luminosities 
j, L 24 M m and using the Bavouzet et al. (2008) vL^^- 

^boi. relation: 

= 6856 x {vLf^f^ (±54%). (3) 

This relation has been calibrated on a sample of 
Spitzer/MIPS galaxies detected at 24, 70 and 160 fim 
with spectroscopic redshifts z < 0.6. The k-corrections 
of the rest-frame 24 fim luminosities used to calibrate this 
formula depend on galaxy IR SED models. However, 
in contrast to other typical recipes to do the vl?^ ^ m - 
Lj^j conversion, the bolometric IR luminosities used in 
the Bavouzet et al. (2008) calibration have been directly 
calculated from the observed fluxes at 24, 70 and 160 
/im. 

We obtained the UV luminosity at a reference rest- 
frame wavelength A = 2000 A for our galaxies using the 
GALEX near-UV (NUV) data for the COSMOS field 
(Schiminovich et al. 2007), which correspond to an ob- 
served effective wavelength A c ff. = 2310 A. We used the 
Bruzual & Chariot (Bruzual 2007) models to compute 
the corresponding k-corrections. We note that, for most 
of our galaxies, the derived UV SFR are a minor fraction 
of the total SFR (< 10% for ~85% of our star-forming 
galaxies and < 20% for ~95% of them). 

Several emission lines are commonly used in the lit- 
erature to trace star formation activity, even though all 
of them are known to present different caveats. B aimer 
lines are affected by underlying stellar absorption and 
dust extinction. [Oil] A3727and [OIII] A5007are sensi- 
tive not only to dust extinction, but also to chemi- 
cal abundance and the degree of ionisation in the star- 
forming galaxy (see e.g. Moustakas et al. 2006). To 
avoid dealing with the additional uncertainties associ- 
ated with the SFR derived from ionisation lines, we only 
study here the Ha A6563 and H/3 A4861 -based SFR, for 
a comparison with the total (IR+UV) SFR obtained for 
the 24 /mi galaxies. 

We computed the SFR derived from the Ha A6563 line 
for galaxies with zCOSMOS redshifts z < 0.3 with the 
Kennicutt (1998) formula: 

SFR Ba (M yr- 1 ) = 7.94 x 10~ 42 x L Hct , (4) 

where Lh q is the luminosity associated with the flux of 
the HaA65631ine expressed in erg s _1 . This formula is 
also valid for a Salpeter IMF over stellar masses (0.1- 
100) M Q . For galaxies with redshifts 0.2 < z < 0.7, we 
computed the SFR derived from the H/3 A4861 line simply 
using equation (j4]) and assuming a case B recombination 
with temperature T = 10, 000K, i.e an intrinsic ratio 
(HaA6563/H/3A4861)int.=2.87 (Osterbrock 1989). 

The four panels in figure M show the SFR derived 
from Ha A6563 and H/3 A4861 line luminosities compared 
to the UV only and total (IR+UV) SFR. Galaxies identi- 
fied with BLAGN from the optical spectra have been ex- 
cluded from this analysis. We also have excluded 5 other 
sources with Ha A6563 and/or H/3 A4861 measurements 
which are obvious cases of AGN, as determined from 
their high X-ray luminosities Lx > 10 42 ' 5 ergs _1 (see 
Section [5]) . We performed the comparison with the UV 
and (IR+UV) SFR before correcting the spectral lines for 



dust, and only including stellar absorption and aperture 
corrections. This allows us to understand the differences 
in the derived SFR because of dust extinction and to as- 
sess how reliable the SFR obtained from the extincted 
Ha A6563 and H/3 A4861 lines are. In Section[6l we derive 
dust extinction corrections for our galaxies by imposing 
that the Ha A6563 or H/3 A4861 SFR are equal to the to- 
tal (IR+UV) SFR. The dust extinction corrections have 
then been applied for all the subsequent analysis made 
in the following sections. 

Several conclusions can be extracted from inspection 
of figure [H Firstly, the general trend observed is that 
the derived SFR are SFR vy < SFR Vmc < SFR (m+VY) . 
As reddening is larger at UV than at optical wavelengths 
(and even larger than at IR wavelengths), this trend ap- 
pears to be a consequence of the differential dust extinc- 
tion. This is supported by the fact that the SFR ob- 
tained from H/3A4861are in better agreement with the 
UV-only SFR than those derived from Ha A6563 (see left- 
hand side panels of figure [8|) . 

For galaxies with i/L 24 ^ 111 < 10 10 L Q at z < 0.3, the ob- 
served Ha A6563 line produces a SFR that is typically a 
factor ~5 below the total (IR+UV) SFR. The differences 
are more dramatic when using the H/3 A4861 line: for 
galaxies with 10 10 < vLf^ m < 10 11 L at 0.2 < z < 0.7, 
the SFR are underestimated by a factor ~ 10. This in- 
dicates that the effects of dust extinction are very im- 
portant in most IR sources. This is consistent with the 
results of previous spectroscopy studies of IR-luminous 
galaxies (e.g. Flores et al. 2004; Liang et al. 2004; Choi 
et al. 2006). The exact factors obtained for the optical- 
to-total SFR conversion mainly depend on the recipes 
adopted to convert mid-IR into bolometric IR luminosi- 
ties. 

For some galaxies in our sample, especially among the 
z/L 24/im > 10 11 L population, the discrepancies between 
the H/3 A4861 -derived and total SFR are much larger 
than the average ( > 50). This could indicate the pres- 
ence of either very-obscured star-forming galaxies or the 
presence of AGN which are to some extent responsible 
for the IR emission. We note that, although we excluded 
from this analysis all the zCOSMOS identified BLAGN, 
narrow- line AGN or mixed AGN/star-forming systems 
might still be present in our sample. We mentioned in 
Section [3] the existence of some IRAC power-law sources 
within our sample. However, only one of the sources with 
highly discrepant SFR is actually an IRAC power-law. In 
some other cases, high ionisation line ratios are observed 
in their optical spectra. As we will discuss in Section [TOl 
star formation and AGN activity might coexist in some 
of these systems. 

6. DUST EXTINCTION 

The discrepancies between the SFR derived from the 
emission lines in the optical spectra and the (IR+UV) 
indicators make necessary the study of dust extinction 
in our galaxies. The common procedure in the literature 
is to adopt a known reddening law to deredden the opti- 
cal spectra. Either the Calzetti et al. (2000), the Milky 
Way (MW; Fitzpatrick 1999 and references therein) or 
the Small Magellanic Cloud (SMC; Prevot et al. 1984) 
reddening laws are generally used. We can test the ade- 
quacy of these different laws to describe the interstellar 
extinction in a sub-sample of our galaxies. 
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We analysed the IR, UV and zCOSMOS spectroscopic 
data of those galaxies in our sample for which both 
Ha A6563 and H/3 A4861 were present and had sufficiently 
reliable measurements, i.e. we restricted the analysis to 
the 24 24 (im galaxies at 0.2 < z < 0.3 with EW> 5 A for 
both Ha A6563 and H/3 A4861 . 

Figure [5] shows the observed Balmer decrement 
L(Ha A6563 ) / X(H/3 A4861 ) of these galaxies, versus the 
logarithm of the ratio between the (far-IR+UV) and 
the UV luminosities. The latter quantity is related to 
the extinction at UV wavelengths ^2000, i-e. (£70^111 + 
^20oo)/^20oo ~ I0°- 4j4200 ° (assuming a uniform geome- 
try for the dust distribution; see e.g. Buat et al. 2002). 
To estimate the far-IR contribution, we used the rest- 
frame 70 /im luminosity, which we compute based on the 
rest- frame 24 /im luminosities and assuming the vl?^ Atm - 
vLl°^ m relation calibrated by Bavouzet et al. (2008). 

The different thick colour lines indicate the relations 
obtained from different reddening laws. We see that both 
the MW and the Calzetti et al. (2000) reddening laws 
produce basically the same curve. This is because of the 
following: we assumed that the observed Ha A6563 -to- 
H/3 A4861 luminosity ratio is given by 



L(Ra) 



2.87 x 10 



-0A(A Ha -A a ) 



(5) 



L(H/3) 

i.e. the intrinsic case B recombination line ratio at- 
tenuated by the corresponding difference of extinctions. 
The extinction ^4h q ,h 3 is proportional to the reddening 
fc(H a ,H^), and thus 

fc(H a )-k(H^) 
A Sa -A^=A 2000 x fc(2oooA) . (6) 

This means that the attenuation on L(Ha)/L(H/3) is 
given by the slope of each reddening curve between the 
HaA6563-H/3A4861and the 2000 A spectral regions. 

As the MW reddening law is characterised by the pres- 
ence of a graphite bump around 2 175 A that is absent in 
the Calzetti et al. law, it might seem surprising that 
both curves in figure [9] are nearly coincident. This is due 
to the wavelengths we considered to study the attenua- 
tion: the slopes of the two reddening laws are basically 
the same between these precise wavelengths. If we had 
used a different UV wavelength for our analysis, the two 
curves would be slightly different. 

Both the MW and the Calzetti et al. reddening laws 
appear to correctly describe the average extinction ob- 
served in our 24 (im galaxies, within the error bars. The 
SMC law, instead, has a quite flatter slope and tends 
to underestimate the observed L(Ha)/L(H/3) ratio. So, 
to deredden the spectral lines in our galaxies, we de- 
cided to adopt the Calzetti et al. reddening law, but 
the MW law should also be suitable for the H/3A4861- 
to-Ha A6563 spectral region (as, in practice, the region 
of the graphite bump is convolved with the NUV filter 
response function, so we are not sensitive to this struc- 
ture). 

We observe in figure [9] that six galaxies in our sample 
display L(Ha)/L(H/3) ratios that are equal to or lower 
than the case B recombination, but at the same time have 
significantly large (£70 + £20oo)/£20oo ratios. Five 
out of six of these galaxies are secure one-to-one 24 fj,m - 
zCOSMOS identifications, so this behaviour is not due to 



incorrect associations. The L(Ha)/L(H/3) ratios of these 
sources can neither be explained by a linear combination 
of two components, one completely free of dust (in a 
case B recombination) and the other one following e.g. 
the Calzetti et al. reddening law (dot-dashed green line 
in figure [9|). 

Instead, the behaviour of the six outlier sources is prob- 
ably the combination of two factors: 1) at 0.2 < z < 0.3, 
the 1-arcsec wide VIMOS slits only cover a small region 
within the galaxies and 2) the dust is not homogeneously 
distributed; in particular, the optical spectra map regions 
where the UV/optical photons can escape with virtually 
no absorption. 

The HST/ ACS /-band images show that all five cases 
with secure identifications are galaxies with a bright op- 
tical nucleus and a quite fainter disky structure around 
(see figure [TU]). This fainter structure is probably ex- 
tracted by dust, which is itself responsible for the IR 
emission. The remaining sixth case corresponds to a 
merger, and the IR emission is likely to be produced by 
the dust surrounding the two merging galaxies. 

Except for these few cases with resolved internal re- 
gions, the application of a single reddening law should 
be suitable for the majority of our galaxies. Although 
sources with non-homogeneous dust distributions might 
also be present at higher redshifts, the possibility of re- 
solving different regions within a single galaxy is much 
less likely. 

Figure [TT] shows the U-band extinctions derived for dif- 
ferent IR luminosity galaxies at different redshifts in our 
sample, assuming the validity of the Calzetti et al. (2000) 
reddening law. To obtain the extinction, we imposed the 
SFR derived from the Ha A6563 or H/3 A4861 line flux to 
be equal to the total (IR+UV) SFR in each galaxy. As 
before, zCOSMOS BLAGN are excluded from this anal- 
ysis. 

We observe an evolution of the median extinction val- 
ues as a function of IR luminosity and redshift. The 
vLl 4 ^ m < 10 10 L Q galaxies within our sample have a me- 
dian extinction Ay — (1.98±0.73) mag, at a median red- 
shift z = 0.22. Among the vhf^ > 10 11 L Q galaxies, 
the median extinction is Ay = (2.70 ± 0.75) mag, at a 
median redshift z = 0.67. These values are consistent 
with previous spectroscopic studies of IR galaxies (e.g. 
Liang et al. 2004; Choi et al. 2006). 

For a comparison, we included in figure [TT1 the extinc- 
tions obtained by Maier et al. (2005) from the spectra 
of 30 Canada-France Redshift Survey (CFRS; Lilly et 
al. 1995) sources at 0.5 < z < 0.9. We see that most 
of the CFRS sources have smaller extinctions than our 
galaxies at similar redshifts. This is expected, as only a 
few of these 30 CFRS sources are IR luminous (Flores et 
al. 1999). 

Choi et al. found that the extinction values within 
their IR galaxy sample were mainly correlated with IR 
luminosity rather than redshift. To test this property 
on our 24 /im galaxies, we computed the correlation co- 
efficients between extinction and IR luminosity, and be- 
tween extinction and redshift: 

^ q )4f W '-^'-' 1 . (7) 



v u y 



where the variable y alternatively indicates rest-frame 
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24 /im luminosity vL v Mm or the redshift z. The symbols 
with an overline indicate mean values and the different a 
correspond to one- variable dispersions. The sum is made 
over all the galaxies considered to probe the correlations. 

We firstly restrict this measurement to a given lumi- 
nosity bin (10 10 < vLl A ^ m < 10 n L Q ) and study the 
correlation between Ay and z. We find that the correla- 
tion coefficient is p(Ay;z) m —0.08 (for N — 67 galax- 
ies), i.e. there is basically no correlation between the two 
variables Ay and z. 

Instead, we can fix the redshift bin and study the cor- 
relation between Ay and vl?^ A ' m . For the 93 galax- 
ies at 0.2 < z < 0.45, we find p{A v ;vLf^ m ) « 0.07. 
For the 48 galaxies at 0.45 < z < 0.70, we have 
p(A v ;vLl ilim ) w 0.46. Thus, within our sample, A v 
appears to be somewhat correlated with IR luminosity 
for the intermediate redshift 0.45 < z < 0.70 galaxies (at 
redshifts comparable to those explored by Choi et al.). 

7. LINE DIAGNOSTICS AND METALLICITIES 
7.1. [OIII] A5007/ H/3 A4861 

We further explored the optical spectral characteris- 
tics of the 24 pm galaxies in COSMOS by analysing the 
variation of typical line ratios as a function of IR lumi- 
nosity and redshift. Whenever possible, we used these 
line ratios to estimate galaxy metallicities. However, in 
contrast to metal abundances, direct line ratios have the 
advantage of being quite model-independent and free of 
physical assumptions. 

Given the wavelength range covered by the zCOSMOS 
spectra, H/3A4861and [OIII] A5007can be measured for 
galaxies at 0.2 < z < 0.7 without being significantly 
affected by fringing. The left panel of figure [12] shows 
the [OIII] A5007 /H/3 A4861 ratios as a function of red- 
shift. Filled circles correspond to our measurements on 
the individual spectra, in the cases when the two line 
EW were > 5 A. The size and colour code for the circles 
is the same as in previous plots. Small empty circles in- 
dicate tentative [OIII] A5007 /H/3 A4861 measurements on 
spectra with H/3A4861EW > 5 A but [OIII] A5007 EW 
< 5 A. The zCOSMOS BLAGN are excluded from this 
analysis. A single case of the shown line ratios corre- 
sponds to an obvious narrow-line AGN, given its high 
X-ray luminosity Lx > 10 42 5 ergs _1 (filled circle with 
a plus sign within) . All line measurements include aper- 
ture and extinction corrections, although these correc- 
tions have little impact on the derived ratios. H/3 A4861 is 
also corrected for stellar absorption in all cases. All the 
completely corrected flux measurements in the cases that 
both lines have EW > 5 A are listed in table [2] 

To study the average line ratios of all our galaxies in- 
dependently of the selection effects imposed by the 5 A 
EW cut, we measured the line fluxes on average compos- 
ite spectra in different IR luminosity and redshift bins 
(asterisks in figure [T2")) . To construct each composite, we 
considered all the galaxies in each bin, independently of 
the line EW measured on their individual spectra. We 
assigned to each composite a redshift equal to the median 
of the individual galaxy redshifts. 

Also, for a comparison, we show the 
[OIII] A5007 /H/3 A4861 ratios of all the zCOSMOS 
galaxies with H/3A4861EW > 5 A (Lamareille et 
al., in preparation; brown dots in figure I12p . The 



average of these ratios are remarkably constant with 
redshift, with log 10 ([OIII] A5007 /H/3 A4861 ) « -0.5 at 
z = 0.2 - 0.3 and log 10 ([OIII] A5007 /H/3 A4861 ) « -0.4 
at z = 0.6 - 0.7. 

By comparison of the [OIII] A5007 /H/3 A4861 line ra- 
tios among our galaxies with different luminosities and 
redshifts, and with the line ratios of the entire zCOSMOS 
galaxy population, we can conclude the following: 

• the average [OIII] A5007 /H/3 A4861 values 
obtained from the stacked spectra for 

io 10 L© galaxies at z ~ 0.25 are compa- 
rable to those of 10 10 < vl?^ ^ m < 10 11 L Q sources 
at z ~ 0.4 — 0.7. This indicates that the hardness 
of the ionising flux is comparable in all these 
galaxies. This fact suggests that the fraction of 
massive stars among all the stars being created in 
these galaxies are similar. 

• For the less luminous IR galaxies in our sample 
at 0.30 < z < 0.45, the average ionisation level 
is quite low and comparable to the average of the 
entire zCOSMOS sample. We obtained the aver- 
age of our IR galaxies by stacking 38 galaxies with 
vjffiB*. < io 10 L Q at redshifts 0.30 < z < 0.45. 
This fact indicates that the ionising fluxes of these 
sources of moderate IR luminosity are similar to 
those of any other optically bright galaxy at low 
redshift. 

• the average [OIII] A5007 /H/3 A4861 ratio obtained 
from the composite spectrum of 42 galaxies with 
i/L^ m > 10 11 Lq at a median redshift of z = 0.67 
is more than 0.6 dex higher than the average ra- 
tio of all the zCOSMOS galaxies at similar red- 
shifts. This result shows that the most lumi- 
nous IR galaxies at z = 0.6 — 0.7 have differ- 
ent chemical properties from other more typical 
optically-selected sources. In particular, among 
our vLl 4f * m > 10 11 L© sources at z > 0.45 with 
H/3A4861and [OIII] A5007EW > 5 A , more than 
a half have log 10 ( [OIII] A5007 /H/3 A4861 ) > 0.5. 
These very high ionisation levels -characterising 
only ~1% of the total zCOSMOS galaxies at these 
redshifts- could suggest the presence of some re- 
maining (narrow-line) AGN within our sample, but 
could also indicate that a large fraction of young 
massive stars of types O and B are being created 
in these galaxies. 

We also compared our derived line ratios with other 
measurements available in the literature. Small green 
and large red crosses in the left-hand panel of figure 
HUshow the [OIII] A5007 /H/3 A4861 line ratios obtained 
by Rupke et al. (2008) on a sample of IRAS-selected 
LIRGs and ULIRGs, respectively, mostly at z < 0.3. 
None of the sources in Rupke et al. have the high ioni- 
sation levels characteristic of several of our vL^ A ^ m > 
10 11 L galaxies at 0.5 < z < 0.7 and a few of the 
10 10 < vLl A » m < 10 11 L galaxies at lower redshifts (al- 
though Veilleux et al. 1995 did find a few cases of ULIRGs 
with high ionisation levels in the local Universe). 

The high ionisation levels characterising some of our 
Pl^vm > io 11 L Q sources at 0.5 < z < 0.7 are compa- 
rable to those found in many sub-millimetre galaxies at 
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z > 1.5 (Takata et al. 2006; right-hand panel of figure 
[T2"|) . Unfortunately, the determination of the source of 
ionisation in both these sub-millimetre galaxies and our 
intermediate-redshift ULIRGs requires much better line 
sensitivity than that provided by low-to-medium resolu- 
tion spectra. Future follow up of these sources with near 
or mid-IR spectroscopy should help us to better constrain 
the nature of our most luminous IR galaxies at interme- 
diate redshifts. 

7.2. [Nil] A6584 /Ha A6563 and the BPT diagram 

For galaxies up to redshift z ~ 0.3, the 
Ha A6563 and [Nil] A6584 lines are observed in the zCOS- 
MOS spectra and quite unaffected by fringing. The 
completely-corrected Ha A6563 and [Nil] A6584 fluxes for 
our 24 \xm galaxies with EW > 5 A are presented in table 
[3] and the corresponding [Nil] A6584 /Ha A6563 ratios are 
shown in figure 1131 We see that our 24 /im galaxies span 
a wide range of [Nil] A6584 /Ha A6563 ratios. This is in 
part due to the same spatial sampling effect discussed 
in Section |4l which produced the large dispersion in the 
Ha A6563EW. 

For ~ 15% of our 24 /im galaxies at z < 0.3, the 
ratios are quite high log 10 ([NII] A6584 / HaA6563) > 
—0.2. Although we do not have reliable H/3A4861and 
[OIII] A5007 measurements for most of them, the high 
[Nil] A6584 / Ha A6563 ratios alone are sufficient to 
classify these objects either as Seyferts or as their 
low-luminosity equivalents, i.e low-ionisation nuclear 
emission-line regions (LINERs; see figured!]). These high 
[Nil] A6584 / Ha A6563 values are more common among 
more luminous IR galaxies at similar redshifts (Rupke et 
al. 2008; crosses in figure fT3]). 

Our 24 fim galaxies at 0.2 < z < 0.3 with measure- 
ments for the four H/3 A4861 , [OIII] A5007 , Ha A6563 and 
[Nil] A6584 lines can be located in the Baldwin, Phillips 
& Terlevich (1981; BPT) diagram (figure EJ). Once 
more, we compare our measurements with the line ra- 
tios of other zCOSMOS galaxies at similar redshifts 
(Lamareille et al., in preparation). Nearly all of our 
24 /im galaxies at these redshifts have mid-IR luminosi- 
ties vLl Atim < 10 10 L Q and occupy quite different posi- 
tions in the BPT diagram, similarly to other optically- 
selected galaxies that are below our detection limits in 
the IR. This indicates the quite heterogeneous charac- 
ter of the IR normal galaxy population: any galaxy with 
modest levels of SF ( < 10 — 15M Q yr _1 ) belongs to this 
category, independently of its degree of chemical evolu- 
tion and already-assembled stellar mass. 

The asterisks in figure [T4l indicate our lines ratios mea- 
sured on the average stacked spectra at median red- 
shifts z — 0.25 and ~ 0.35. Although we did not mea- 
sure Ha A6563 and [Nil] A6584 line fluxes at z > 0.3 on 
an individual basis to avoid being severely affected by 
fringing, the signal-to-noise ratio in the stacked spec- 
tra is sufficiently good as to allow for a reliable mea- 
surement of these lines (see figure [T5j ■ These average 
spectra correspond to 10 10 < vLl 4 ^ m < 10 11 L and 
yj^ifim < io 10 L galaxies and their line ratios are lo- 
cated just beyond the empirical line proposed by Kauff- 
mann et al. (2003a), which delimits the transition region 
between pure star- forming and mixed AGN/star- forming 
galaxies in the BPT diagram. This fact suggests that sev- 



eral mid-IR-selected galaxies at z ~ 0.3 — 0.4 redshifts, 
even with modest IR luminosities, might have a compos- 
ite star-forming/ AGN nature. 

7.3. Oxygen abundances 

We used the measured and corrected line fluxes to 
estimate the oxygen abundances [O/H] of our 24 /jm- 
selected galaxies, whenever possible. 

Several methods exist to derive abundances from 
emission-line ratios, all of which have both advan- 
tages and disadvantages. For example, the Pet- 
tini & Pagcl (2004) formulation uses both the 
[OIII] A5007/H/3A4861 and [Nil] A6584 /Ha A6563 ratios, 
but the [Nil] A6584 /Ha A6563 ratio is particularly af- 
fected by AGN radiation fields, so the resulting abun- 
dance can be very uncertain when part of the spectral 
light has a non-stellar origin. The more widely used R23 
parameter (Pagel et al. 1979) is based on the [Oil] A3727 , 
H/3 A4861 and [OIII] AA4959, 5007 lines, all of which are 
usually observed in optical spectra. However, the oxygen 
abundance is two- valued with respect to i?23, and any 
given value of R23 lower than the maximum is compat- 
ible with a low- and a high-metallicity solutions. When 
available, additional emission lines are sometimes used 
to try to lift this degeneracy. Another common approach 
is to consider only one of the branches (usually that of 
high-metallicity) in the [O/H] versus R23 relation. Fur- 
ther discussion about different abundance estimators can 
be found in e.g. Kobulnicky et al. (1999) and Kewley & 
Dopita (2002). 

We used the algorithm described by Maier et al. (2005) 
to derive oxygen abundances, which is based on the Kew- 
ley & Dopita (2002) model. This algorithm also uses the 
information on different line fluxes, but performs a si- 
multaneous fit to all available emission lines, with three 
free parameters: the extinction, the ionisation parame- 
ter (which is proportional to the flux of ionising photons) 
and the abundance [O/H]. 

We implemented this algorithm in two redshift ranges: 

1) 0.2 < z < 0.3, for the same galaxies shown 
in the BPT diagram, for which we have H/3 A4861 , 
[OIII] A5007 , Ha A6563 and [Nil] A6584 measurements; 

2) 0.5 < z < 0.7, for galaxies with [OII]A3727, 
H/3A4861and [OIII] A5007 measurements. In all cases, 
we fixed the value for the extinction to that derived from 
the balance between the Ha A6563 or H/3 A4861 and the 
total (IR+UV) SFR, i.e. we used the algorithm with 
only two free parameters. We excluded from this analy- 
sis a few sources belonging to the group of galaxies with 
H/3 A4861 /Ha A6563 ratios close to a case B recombina- 
tion, but large (£70^111 + ^20oo)/^20oo ratios (see Section 
E). All zCOSMOS BLAGN are also excluded (as they 
were from all our previous line diagnostic). 

For sources in the redshift range 1), the algorithm 
nearly always produces either non-degenerate solutions 
or two solutions of relatively high-abundance, in such a 
way that a useful lower limit can be set. Galaxies in 
range 2), in contrast, have the same problem as when 
using the R23 parameter: unless we are close to a max- 
imum R23 value, there are one possible high- and one 
low-metallicity solutions. 

Figure \W\ shows the obtained abundance values as a 
function of redshift. Non-degenerate cases are indicated 
with filled circles with error bars. Circles with upward- 
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pointing arrows at z < 0.3 indicate lower limits. As 
before, small empty circles represent galaxies for which 
Ha A6563 and H/3 A4861 EW > 5 A, but with at least one 
of the other line EW below that cut. 

Among our 26 galaxies with computed abundances at 
0.5 < z < 0.7, 9 (~35%) have non-degenerate values. 
For the remaining 17, there is an upper-branch (large 
empty circles) and a lower-branch (large black crosses) 
solutions. The dashed line indicates the solar abundance 
derived by Asplund et al. (2004). As we have computed 
all the abundances using the same technique, it makes 
sense to compare the results we obtained at different red- 
shifts. 

The ~ 35% of non-degenerate abundances at redshifts 
0.5 < z < 0.7 are solar or sub-solar, while most of the 
lower redshift metallicities at 0.2 < z < 0.3 are so- 
lar or super-solar. We note that the important point 
here is the relative values of the abundances at differ- 
ent redshifts. The fact they are sub- or super-solar 
is incidental in our comparison, as actually different 
metallicity indicators have systematic differences of up 
to ~ 0.5 dex among them (e.g. Kewley & Dopita 2002; 
Ellison & Kewley 2006). For example, if we used the 
Pettini & Pagel (2004) prescription for our galaxies at 
0.2 < z < 0.3, we would obtain abundances which are 
~ 0.2 — 0.3 dex lower than those shown in figure [TBI 

At 0.5 < z < 0.7, non-degenerate abundances 
correspond to cases with large R23 values, i.e. 
with relatively high [OIII] A5007 /H/3 A4861 and 
[Oil] A3727 /H/3 A4861 ratios. These ratios are illus- 
trated in figure [T7] We see that all these galaxies with 
non-degenerate abundances lie within or very close 
to the star-forming galaxies/ AGN separation region 
proposed by Lamareille et al. (2004) . This is a potential 
concern in our analysis: the derived metallicities of 
these objects could be affected by the plausible presence 
of an AGN. We note, however, that there is only 
one X-ray AGN among these galaxies and it lies just 
below the separation region (red large circle behind the 
green small one), showing that the association of high 
ionisation lines and AGN activity is not straightfor- 
ward and has to be taken with care. Besides, strong 
ionisation lines such as [Oil] A3727 and [OIIIJA5007- 
resulting in very high [OIII] A5007 /H/3 A4861 and 
[Oil] A3727 /H/3 A4861 values- can also be produced 
by very young massive stars (see e.g. Chariot & 
Longhetti 2001), with no need of invoking the presence 
of an AGN. 

The galaxies with degenerate abundances at 0.5 < 
z < 0.7 are likely to have high metallicities (see the 
mass-metallicity relation discussion in Section 17.41 and 
figure [T5|) . They constitute ~65% of our higher-redshift 
metallicity sample. These upper-branch metallicity val- 
ues are comparable to those we obtained at lower-redshift 
sources and those derived for JiM.S'-selected LIRGs and 
ULIRGs by Rupke et al. (2008; small green and red 
crosses in figure [T6|) . The remaining 35% of our galaxies 
at 0.5 < z < 0.7 have lower abundances than the vast 
majority of low-redshift IR sources. These results sug- 
gest that at least 1 /3 of the hosts of the most intense IR 
emission have chemically evolved from redshifts z ~ 0.7 
toz~ 0.2. 

7.4. The mass-metallicity relation 



The mass (or luminosity)-metallicity relation for differ- 
ent galaxy populations has been the subject of multiple 
studies in the literature (e.g. Zaritsky et al. 1994; Con- 
tini et al. 2002; Lilly et al. 2003; Tremonti et al. 2004), 
and it has been found to evolve with redshift (e.g. Kob- 
ulnicky et al. 2003; Contini et al. 2008). This relation 
can be explained within the context of hierarchical mod- 
els of galaxy formation, and by including the effects of 
supernovae-driven winds which help to remove metals es- 
pecially out of the less-massive galaxies (e.g. Tissera et 
al. 2005; Brooks et al. 2007). 

It is then of interest to see whether IR-selected galax- 
ies follow a similar relation. Figure [18] shows metal- 
licities versus already-assembled stellar masses for our 
24 /im galaxies (at 0.2 < z < 0.3 and 0.5 < z < 0.7, in 
the top and bottom panels, respectively). Stellar masses 
for all zCOSMOS galaxies have been computed by fit- 
ting Bruzual & Chariot models to the multiwavelength 
(U-band to 4.5 jum) SED of these sources (Bolzonella et 
al., in preparation). All stellar masses assume a Salpeter 
(1955) initial mass function (IMF) over stellar masses 
M = (0.1 — 100) Mq. As before, filled circles correspond 
to non-degenerate metallicity values, while empty circles 
and crosses show the upper and low-branch solutions, 
respectively, in all the degenerate cases. For a compar- 
ison, we have also included a sample of 0.5 < z < 0.9 
galaxies from the CFRS (triangles in the bottom panel 
of figure [18]) . All these galaxies have secure metallicity 
determinations (Maier et al. 2005). Three of them have 
been detected by ISO at 15 /im (Flores et al. 1999; yellow 
filled triangles). 

The solid line in both panels of figure [18] is the best 
mass-metallicity relation fitted by Tremonti et al. (2004) 
on Sloan Digital Sky Survey (SDSS) galaxies at z ~ 0.1. 
The dashed lines in the bottom panel are the correspond- 
ing relations at 0.5 < z < 0.7 and 0.7 < z < 1.0, as 
obtained by Contini et al. (2008) based on the entire 
zCOSMOS-bright 10k sample. Both the Tremonti et al. 
and the Contini et al. relations have been re-scaled to 
the Salpeter IMF. Instead, we did not modify any of their 
metallicities. The Tremonti et al. calibration produces 
very similar metallicities as the Kewley & Dopita method 
(Ellison & Kewley 2006). The Contini et al. metallicities 
are basically the same as ours for the common zCOSMOS 
objects. Thus, we left the Contini et al. metallicities in 
their original scale. 

Inspection of the top panel of figure [18] shows that 
some of our 24 /im -selected galaxies at 0.2 < z < 0.3 
have metallicities which are at most marginally consis- 
tent with the general mass-metallicity relation at z ~ 0.1. 
However, in many cases, we only show lower limits, so 
they still have the possibility of being in agreement with 
the general mass-metallicity trend. 

At higher redshifts, ^65% of our galaxies have one of 
the two possible metallicity values consistent with the 
general mass-metallicity relation derived by Contini et 
al. (2008) for all zCOSMOS-bright 10k galaxies. How- 
ever, most of our galaxies with a non-degenerate metal- 
licity determination appear to have lower abundances 
than expected. Only 2 out of these 9 secure-metallicity 
galaxies have stellar masses and abundances more or less 
in consistence with the extrapolation of the Contini et 
al. relations. The other 7 (i.e. ~ 27% of the total 
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0.5 < z < 0.7 sample) are significantly out of the general 
mass-metallicity trend. 

Several of the CFRS galaxies shown in figure [15] have 
also been found to be under-abundant with respect 
to galaxies in the upper metallicity branch (Maier et 
al. 2005). We note, however, that their characteristic 
stellar masses are, in general, quite smaller than for most 
of our galaxies. So, the discrepancies we observe for some 
of our IR galaxies are considerably more important. 

As we discussed in Section 17.31 most of these galax- 
ies have high ionisation levels, so there exist the concern 
that the plausible presence of a narrow-line AGN could 
be in part responsible for the low measured metallicitics. 
However, it is not clear that the discrepancies with the 
general mass-metallicity relation are due to the effects 
of a plausible AGN. For example, one of the sources 
with the highest ionisation ratios logio(OIII/H/3) = 
0.90 ± 0.10 (see figure [T7|) has a stellar mass consis- 
tent with the extrapolation of the Contini et al. re- 
lations. On the other hand, two out of three of the 
galaxies with secure metallicities that are just within 
the star burst region in the [OIII] A5007 /H/3 A4861 versus 
[Oil] A3727 /H/3 A4861 diagram have metallicities signifi- 
cantly below the expected trend. 

From the study of a sample of IRAS-selected galax- 
ies, Rupke et al. (2008) found that the vast majority of 
LIRGs and ULIRGs at low z < 0.3 redshifts are signif- 
icantly under- abundant. We cannot confirm this from 
our metallicity measurements at similar redshifts, but 
the difference is probably due to the fact that our IR 
galaxies at z < 0.3 are less IR-luminous than Rupke et 
al. galaxies and more similar to typical optical galaxies 
at low redshifts. 

Instead, we can confirm the under-abundance for a 
fraction (~ 1/3) of our 24 /urn sources at 0.5 < z < 0.7. 
To explain the relatively low metallicities of their galax- 
ies, Rupke et al. proposed that, during the IR phase, 
less-abundant gas from outer parts of the galaxies could 
flow inward to the central regions and dilute the nuclear 
metal content. Although this mechanism appears to be 
common in gas-rich mergers (e.g. Barnes & Hcrnquist 
1996), its influence is less clear in regular star- forming 
disks, as many luminous IR galaxies at z < 1 are known 
to be. Supernovae-driven outflows are also known to 
effectively remove metals from a galaxy (e.g. Dalcan- 
ton 2007) but it seems implausible that this mechanism 
alone can explain the under-abundance observed in IR 
galaxies (see the discussion in Rupke et al. 2008). More 
likely, a combination of the two phenomena could have 
taken place at the same time. Thus, gas mixing produced 
by inflows and outflows during the IR phase could tem- 
porarily move galaxies appart from the mass-metallicity 
relation. 

8. AGN IN THE 24 fim SAMPLE 

Optical spectra also allow us to test the presence of 
AGN within our 24 /im sample. The zCOSMOS clas- 
sification includes a special flag for BLAGN, gener- 
ally recognised by the presence of broad Ha A6563 or 
H/3 A4861 lines, for the shape of the continuum, and also 
for the presence of Mgll and CIII in emission in high-z 
sources. 

The entire zCOSMOS-bright 10k catalogue contains 
132 BLAGN with good-quality spectroscopic redshifts 



and spectra, all of which have I < 22.5 AB mag. 
64 out of these 132 BLAGN are 24 /mi -detected with 
5 , 24^ m > 0.30 mjy. The histograms in the left panel of 
figure [19] show the redshift distributions of the two sam- 
ples, both of which include objects up to z > 3. 

The comparison of both redshift distributions shows 
that, at least up to z « 2, the fraction of I < 22.5 
BLAGN that are IR-detected has little variation with 
redshift. This fact indicates that the sensitivity limit 
of the 24 jjLUY survey is not the only factor that deter- 
mines the IR detection of these sources. Other physical 
differences might mean that some BLAGN are not IR- 
detected, even when lying at lower redshifts than others 
that are detected. 

Within the unification scheme for AGN (An- 
tonucci 1993), broad lines are present in the spectra when 
the source is oriented face-on and the central gas clouds 
are visible. The surrounding dusty torus is responsible 
for the IR emission. The flux at IR wavelengths depends 
on the incoming flux of UV photons from the central en- 
gine, so it is related to the amount of material accreted 
onto the central black hole. Thus, the difference in the 
IR emission of these type-1 AGN could be directly indi- 
cating the accreting power in these objects. 

As expected, the most luminous IR sources in our flux- 
limited 24 /im sample are AGN at redshifts z > 1 (cf. 
figure [3]), particularly BLAGN. AGN are known to dom- 
inate the bright-end of the mid-IR luminosity function 
at high redshifts (Caputi et al. 2007). To investigate 
the presence of other obvious (narrow-line) AGN within 
our sample, we looked for X-ray counterparts with Lx > 
10 42.5 ergs -i i n XMM-Newton maps of the COSMOS 
field (Hasinger et al. 2007). We found 35 IR sources as- 
sociated with narrow-line Lx > 10 42 5 ergs _1 sources. 
The distribution of 24 /im luminosities vl?^ Mm of these 
galaxies compared to the distribution of luminosities for 
the BLAGN is shown in the right panel of figure [T9] It 
is clear from these two distributions that the BLAGN 
in our sample are considerably more luminous in the IR 
than the most powerful narrow-line AGN. 

The fact that all of the identified most-luminous IR 
AGN are BLAGN might be partially due to the opti- 
cal magnitude cut limiting the zCOSMOS-bright survey. 
Some narrow-line AGN at high redshifts could also be 
IR-luminous, but they are probably too faint in the op- 
tical bands to be zCOSMOS targets. The prevalence of 
BLAGN among the most-luminous IR galaxies detected 
in zCOSMOS-bright could also be consequence of the ori- 
entation effect: in BLAGN, the optical nucleus is more 
exposed than in edge-on-oriented active sources. 

On the other hand, 45 of our 24 /im sources are charac- 
terised by an IRAC power-law SED, as defined in Section 
[3] 30 out of 45 of these sources are also zCOSMOS- 
bright BLAGN. To inspect whether the remaining 15 
IRAC power-law sources had also an independent sig- 
nature of AGN activity, we looked for counterparts of 
these sources in the X-ray catalogues for the COSMOS 
field (Hasinger et al. 2007). We found that 11 out of 15 of 
the IRAC power-law sources not identified as zCOSMOS 
BLAGN are detected in X-rays. 

The composite zCOSMOS spectra of these objects 
also suggests the presence of active sources. Figure [20] 
shows the composite average spectrum of 8 IRAC power- 
law sources at 0.5 < z < 0.85 for which the individ- 
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ual zCOSMOS spectra do not correspond to BLAGN. 
We see in this spectrum high [OIII] A5007 /H/3 A4861 and 
[OIII] A5007 / [Oil] A3727 line ratios and also the presence 
of [Nelll] in emission. All these characteristics suggest 
that most of the IRAC power-law sources which are not 
individually identified with zCOSMOS BLAGN are ei- 
ther narrow-line AGN or composite AGN/starburst sys- 
tems. 

9. THE INFRARED PHASE WITHIN THE GALAXY 
STAR-FORMATION HISTORY 

Optical spectra have also features that provide in- 
formation on the past star-formation history (SFH) of 
galaxies. For example, the presence of higher-order 
Balmer lines in absorption (see figures 0] to [S]) indi- 
cates recent activity, i.e. star formation occurred within 
the last Gyr (e.g. Dressier & Gunn 1983). Recently, 
Kauffmann et al. (2003b) studied 4000 A-break strengths 
D n (4000) and H<5 EW for a sample of local SDSS galaxies 
and showed how these two quantities can be used to put 
constraints on the galaxy mean stellar age and SFH. 

Figure [2D shows the expected evolution of the H<5 EW 
versus D n (4000) as a function of age for galaxies with 
different SFHs (in this section, positive values of the EW 
indicate line absorption). We modelled this evolution 
using the Bruzual & Chariot (Bruzual 2007) templates 
for galaxies with exponentially-declining star formations 
(with r = 0.01 and 0.1 Gyr; solid and dotted lines, re- 
spectively) and a constant star formation (dashed line). 
We measured D n (4000) as the ratio of the average flux 
densities between the (4000 - 4100) and (3850 - 3950) A 
bands (Balogh et al. 1999). We only considered mod- 
els with solar metallicities, as both D n (4000) and the 
Hi5 EW have little dependence on metallicity at ages 
< 1 Gyr after the onset of star formation (see Kauffmann 
et al. 2003b). As we show below, these are the timescales 
of relevance for IR galaxies. 

Our aim is then to locate IR galaxies in the US EW - 
D n (4000) diagram and see whether we can constrain at 
what stage of a galaxy life the IR phase is occurring. 

The measurement of D n (4000) can be easily done on 
individual spectra. We show in figure [22] the result- 
ing D n (4000) measurements for all our 24 /im galaxies 
at 0.6 < z < 1.0. The zCOSMOS BLAGN have been 
excluded from this analysis. We corrected the D n (4000) 
measurements for dust extinction, although these correc- 
tions are very small, given the narrow wavelength range 
considered. For a comparison, we show the D n (4000) 
measurements for all the non-IR-detected zCOSMOS 
galaxies in the same redshift range (brown dots in fig- 
ure!^ Lamareille et al., in preparation). 

We see from figure H21 that the D n (4000) strengths 
of IR-selected galaxies are in all cases quite small 
(D n (4000) < 1.3), and they do not depend on redshift. 
The values we obtain for D n (4000) in 24 /im galaxies 
are in agreement with those measured by Marcillac et 
al. (2006) on a sample of 25 /S'O-selected LIRGs at 
z ~ 0.7. The small values for D n (4000) indicate the 
presence of young generations of stars in IR galaxies. We 
note that IR galaxies display the smallest D n (4000) val- 
ues among zCOSMOS galaxies, for which D n (4000) can 
take values up to 2 or more. 

In fact, the existence of young stars is expected in 
galaxies that are IR-selected, i.e. that are selected to 



have on-going star formation. However, an old (> 1 Gyr) 
galaxy with on-going declining star formation displays 
quite redder colours. The fact that D n (4000) is invari- 
ably small indicates that IR galaxies either are young 
systems or are experiencing a rejuvenation of their stel- 
lar content. 

Unfortunately, the US EW are quite small and the as- 
sociated errors too large as to rely on individual mea- 
surements in most cases. We then measured US only on 
average stacked spectra, for which the signal-to-noise ra- 
tio is much higher. We constructed 11 composites corre- 
sponding to galaxies of similar IR luminosities in narrow 
redshift bins between z = 0.6 and 1.0. Each composite 
is the average of 17 to 32 individual spectra, depending 
on the case. 

To correct the H<5 EW for line filling, we measured 
the H7 emission line on each composite, corrected it for 
stellar absorption and assumed an intrinsic decrement 
H7/H15 =1.80 (for a case B recombination with temper- 
ature T = 10,000K; Osterbrock 1989). Beyond z ~ 0.9, 
H7 is within the fringing-affected region of our spectra 
and its measurement becomes quite uncertain. Thus, at 
these high redshifts, we simply used the same line-filling 
corrections as in immediately lower redshift bins. The 
resulting average US EW versus D n (4000) measurements 
are shown with asterisks in figure 1231 

IR galaxies occupy, on average, a well-defined place 
in the H<5 EW - D n (4000) diagram. The ideal single 
constant SFH is consistent with the average location of 
some IR galaxies but, still in this case, the associated 
ages are very young (< lGyr). For most IR galaxies, 
however, the track of a constant SFH is well above their 
location in the US EW - D n (4000) diagram. 

The locus of IR galaxies in the RS EW - D n (4000) plane 
can, instead, be well explained by secondary bursts of 
star formation (thin solid lines in figure[23j see also Ham- 
mer et al. 2005; Marcillac et al. 2006). These secondary 
bursts only need to form an additional minor amount 
(5-10%) of stellar mass in these galaxies to produce the 
Hi5 EW and D n (4000) values observed. Here we confirm 
that such a scenario can explain the average properties of 
the brightest IR galaxies found all the way from redshifts 
z — 0.6 through 1.0. 

Two additional interesting conclusions can be ex- 
tracted from figure l2lfl firstly, the small D n (4000) values 
and relatively small US EW displayed by some compos- 
ites indicates that the rejuvenation of the stellar popu- 
lations occurs, in many cases, in relatively old (> 1 Gyr) 
galaxies. We know that, at z < 1, LIRGs and ULIRGs 
are generally associated with intermediate-to-large stel- 
lar mass galaxies (M > 1 x 10 10 M Q ; see Caputi et 
al. 2006a). So, the brightest IR emission at these red- 
shifts is being hosted by galaxies which already formed 
most of their stellar mass at earlier epochs. 

Secondly, the average location of our sources on the 
US EW - D n (4000) plane strongly constrains when the 
IR phase is being produced: on average, galaxies are IR 
bright in the elapsed time between 10 7 and 10 8 years 
after the onset of the secondary burst. We note that 
our results do not completely exclude the possibility 
that some galaxies have more prolonged IR phases char- 
acterised by more or less constant SFR. Galaxies that 
still have sufficiently large gas reservoirs remaining at 
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z < 1 could in principle sustain the typical SFR of LIRGs 
(~ 10 — 100Mo/yr) for some longer times. However, 
given the already-assembled stellar masses of LIRGs and 
ULIRGs, this situation is probably implausible in most 
cases. 

Most of the dust contained in IR galaxies is produced 
by AGB stars, which start to inject dust into the in- 
terstellar medium only a few hundreds of Myr after the 
production of the burst (e.g. Dwek 2005). The fact that 
the IR phase occurs between 10 7 and 10 8 years after the 
onset of the secondary burst indicates that the dust be- 
ing heated during the IR phase is probably 'old dust', 
i.e. dust that has been created by previous generations 
of stars. 

10. YOUNG POWERFUL STARBURSTS OR FURTHER AGN 
ACTIVITY? 

We have seen in the previous section that the bulk 
of the IR phase in galaxies at z < 1 is observed, on 
average, 10 7 — 10 s years after the onset of a secondary 
burst of star formation. The minimum average time for 
the detection of the IR phase is in part simply due to 
the lower rate of galaxies that can be detected in shorter 
(< 10 7 yr) elapsed times. But, if this were the only factor 
determining the rate of detection of younger starbursts, 
one would possibly expect to find around one out of ten 
IR galaxies being at the earliest stages of a new burst of 
star formation. 

We looked for galaxies with signatures of young (< 
10 7 yr) starbursts in their optical spectra within our 
24 /im sample. At early stages of star formation, the EW 
of absorption features associated with high-order Balmer 
lines is expected to be small (cf. figure [23)) . In addition, 
star formation should be sufficiently strong as to clearly 
detect even the higher-order Balmer lines in emission (ex- 
cept, maybe, for extremely-obscured sources). As we ex- 
plained in Section (9) the typical EW widths measured 
for US on different spectra are, in general, small and the 
associated errors quite large as to consider these measure- 
ments on an individual basis. In spite of this, we used 
these tentative measurements of US to look for candidate 
galaxies to be hosts of very young starbursts. 

We found four likely candidates for young starbursts 
within our 24 /im sample. The spectra of all of these 
galaxies are characterised by the presence of US and He in 
emission. One case also shows a slightly broad Hell line 
(see example in figure [24)) . which is a signature of stellar 
winds originated in Wolf-Rayet stars in young-starburst 
galaxies (Schaerer et al. 1999). 

We note, however, that other characteristics of these 
spectra could also suggest the possibility of AGN ac- 
tivity: in particular, the presence of a [Nelll] emission 
line or the high [OIII] A5007 / [Oil] A3727 ratios. Actu- 
ally, both star formation and AGN activity are likely 
to coexist in these objects. Unfortunately, the lack of 
spectral coverage for the Ha A6563 and [Nil] A6584 lines 
prevents us from further constraining the nature of these 
sources. 

We looked for X-ray counterparts of these objects 
in the XMM-Newton catalogue for the COSMOS field 
(Hasinger et al. 2007). None of them appear to be iden- 
tified with an X-ray source. These galaxies are not iden- 
tified either in the deeper Chandra Observatory X-ray 
maps for the COSMOS field (C-COSMOS; P.I. M. Elvis). 



We found three more galaxies in our sample that also 
have H<5 and He in emission, but the lines are weaker and 
the noise in the spectra is higher. Some extra candidates 
are probably missing because of extinction effects. Still, 
these presumably young starbursts seem to be quite less 
than ~10% of our galaxy sample. This fact suggests that 
the onset of the IR phase within ~ 10 7 yr after the burst 
might be indeed a rare phenomenon. Of course, given the 
optical magnitude cut of the zCOSMOS survey, we can- 
not exclude that IR galaxies at the very earliest stages 
of a rejuvenating starburst are more common among op- 
tically faint sources. 

The low fraction of very young starbursts among IR- 
selected galaxies is a known phenomenon in the local 
Universe. Roussel et al. (2006) studied the case of NGC 
1377, an extremely young ( < 1 Myr) local starburst, se- 
lected for being an outlier in the far-IR/radio correlation. 
These authors note that such galaxies represent ~ 1% of 
faint IRAS galaxies. In contrast to our cases, the opti- 
cal spectrum of NGC 1377 shows very weak [Oil] A3727 , 
[OIII] A5007 and H/3 A4861 emission, due to dust obscu- 
ration. The mid-IR spectrum of this source, instead, is 
dominated by molecular-hydrogen Hi lines, indicating a 
young reservoir of star formation. 

In a forthcoming paper, we will analyse the pres- 
ence of young starburst candidates within the deeper 
24 /im galaxy sample for the COSMOS field up to higher 
rcdshifts. This will allow us to constrain whether a 
fainter IR phase is more often associated with these 
sources at the earliest stages of star formation. 

11. SUMMARY AND CONCLUSIONS 

In this work, we have presented the results of the study 
of zCOSMOS optical spectra for 609 24 /im galaxies with 
"5a4^m > 0.30 mJy, selected over 1.5 deg 2 of the COS- 
MOS field. This is the most extensive analysis of the 
spectra of a very large sample of mid-IR galaxies anal- 
ysed to date that covers the redshift ranges < z < 1 and 
< z < 3 for normal galaxies and AGN, respectively. 

Depending on the lines present in the spectra, we made 
different diagnostics for sources at different redshifts. 
As expected, at low z < 0.3 redshifts, our sample is 
dominated by IR-normal galaxies (in our classification, 
strictly, vLl 4 "™ < 10 10 L o ). The optical spectra of these 
sources display a large variety of emission-line EW and 
ratios. This is in part due to a spatial sampling effect, 
because the 1-arcsec-width VIMOS slit maps only the 
central region within each galaxy, which in some cases 
is not representative of the entire galaxy. The same ef- 
fect has been observed in the study of low-redshift IRAS 
galaxy optical spectra. But the variety of spectral char- 
acteristics also reveals the different physical conditions 
of the hosts of normal IR galaxies. This class includes 
galaxies with different degrees of chemical evolution and 
already-assembled stellar masses: from some low-mass, 
low-metallicity galaxies to many already massive, chem- 
ically enriched, galaxies. 

At higher redshifts, we studied the variation of the ion- 
isation levels in different luminosity IR galaxies, through 
their [OIII] A5007/H^A4861 line ratios. We found that 
the average [OIII] A5007/H/3 A4861 ratio for vLf^ > 
10 11 L/0 galaxies at z ~ 0.6 — 0.7 is more than 0.6 dex 
higher than the average ratio of all zCOSMOS galax- 
ies at similar redshifts. This allows us to conclude 
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that ULIRGs have distinct optical spectral properties 
when compared to typical optically-selected sources. In 
particular, for some of our galaxies, we find very high 
log 10 ([OIII] A5007/H/3A4861) > 0.5 ratios, comparable 
to those of sub-millimetre galaxies at z > 1.5. These 
are cases where massive star formation and narrow-line 
AGN activity are likely to coexist. Near or mid-IR spec- 
troscopy is necessary to provide further information on 
the ionisation source in these galaxies. 

Around 1/3 of the galaxies for which we could derive 
metallicities at 0.5 < z < 0.7 are unambiguously below 
the general mass-metallicity relation at the correspond- 
ing redshift. Although, in some cases, the plausible pres- 
ence of an AGN could affect the derived metallicities, this 
is probably not the main cause of the discrepancies with 
the mass-metallicity relation. Instead, gas mixing due 
to inflows and outflows during the IR phase could tem- 
porarily move sources out of the mass-metallicity trend 
(see e.g. Rupke et al. 2008). 

The comparison of optical emission-line derived and 
fiducial (IR+UV) SFR indicates that IR galaxies are 
generally quite affected by extinction. SFR derived 
from the extincted HaA6563and H/3 A4861 lines under- 
estimate the real SFR, on average, by factors ~ 5 and 
~ 10, respectively. We note, however, that these dis- 
crepancies are expected to be maximal for the most lu- 
minous IR galaxies. In other galaxies, the Ha A6563 or 
H/3 A4861 luminosities should constitute more reliable 
SFR estimators. 

Also from our spectroscopic analysis, we found that it 
is difficult to explain the properties of many IR galaxies 
at z < 1 with a constant SFH. Instead, a SFH charac- 
terised by secondary bursts of star formation can simulta- 
neously explain the blue continuum and relatively small 
HS EW observed in most of these sources (see also Ham- 
mer et al. 2005). These results also constrain the epoch 
of IR activity: the LIRG and ULIRG phases for star- 
forming galaxies occurs, on average, between 10 7 and 10 s 
years after the onset of the secondary burst. Most of the 
hosts of intense IR activity at 0.6 < z < 1.0 have older 
underlying stellar populations and have assembled most 
of their stellar mass in the past (Caputi et al. 2006a). 
The very short times for the onset of the IR phase after 
the burst also suggest that the dust being heated during 
these rejuvenating episodes of star formation must be 
'old dust', i.e. dust produced by AGB stars from previ- 
ous stellar generations. All these conclusions are valid in 
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general for all the galaxies we have analysed here, inde- 
pendently of their IR luminosity and redshift, suggesting 
that the way in which star formation proceeds is basically 
the same in most of these systems (i.e. independently of 
the mechanism triggering the burst). 

Finally, ~10% of our galaxies are recognised as 
BLAGN from their zCOSMOS spectra. Another ~ 
5 — 6% are X-ray-luminous AGN with narrow-line spec- 
tra. These figures are consistent with the fraction of 
AGN-dominated IR sources found in the literature (e.g. 
Alonso-Herrero et al. 2006; Caputi et al.2006a,b,2007). 
However, AGN activity might also be present in some 
other galaxies, as suggested by their high ionisation lev- 
els or IRAC-power law SEDs. The exact fraction of 
AGN and AGN/star- forming composite systems among 
IR galaxies is still far from being clear. In the future, very 
deep X-ray data and far-IR data obtained with the Her- 
schel telescope will be available for the COSMOS field. 
Both datasets analysed in conjunction should help to put 
stronger constraints on the AGN contribution to the IR 
background. 

The results obtained in this paper offer some impor- 
tant clues of the role of the IR phase in star-formation 
history and galaxy evolution, mainly at z < 1 . IR sources 
are the signposts of the most intense star formation ac- 
tivity and on-going build-up of stellar mass. But the IR 
phase is a kind of a transitory state in which some galaxy 
properties (e.g. the strength of the 4000 A break and, in 
some cases, the metallicity) are perturbed from their reg- 
ular values in the IR-inactive galaxy life. This should be 
kept in mind when considering IR sources for a census 
of galaxy properties, as the values measured for some of 
these properties could be well apart from those that will 
hold at long term once the IR phase has faded out. 
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TABLE 1 

THE NUMBERS OF 24 £tm GALAXIES WITH 524 > 0.30 mjy which satisfy the different selection criteria 
IMPOSED BY THE ZCOSMOS-BRIGHT SURVEY. EACH CATEGORY IS INCLUDED IN THE PREVIOUS ONE. ALL 
ASSOCIATIONS HAVE BEEN DONE WITHIN A 2 ARCSEC MATCHING DISTANCE. 



Selection criterion Number of sources 



Sources in SCOSMOS shallow (1.75 x 1.97 deg 2 ) 9,807 

Sources in the zCOSMOS-bright lOk-sample field (1.5 deg 2 ) 3,150 

Sources with / < 22.5 AB mag 2,084 

Sources with associations in zCOSMOS-bright 10k sample 703 

Secure or likely zCOSMOS-bright associations 668 

Sources with secure zCOSMOS-bright spectra 611 
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TABLE 2 

The H/3 A4861 and [OIII] A5007 line fluxes for the galaxies in our sample for which 

BOTH LINE EW ARE GREATER THAN 5 A. ALL FLUXES ARE APERTURE AND 
EXTINCTION-CORRECTED. H/3 A4861 FLUXES ALSO INCLUDE CORRECTIONS FOR STELLAR 

absorption. The flux units are 10 -15 ergs" 1 ™^ 2 . zCOSMOS BLAGN have been 

EXCLUDED. THE SOURCES FLAGGED WITH AN ASTERISK CORRESPOND TO X-RAY AGN WITH 

L x > lO^'ergs- 1 . 



zCOSMOS id 


Zspec 


H/3 A4861 


[OIII] A5007 


\og m uLl 4,lrn 


813537 


0.2602 


1.74±0.42 


1.34±0.41 


< 10 


812917 


0.3398 


1.41±0.28 


1.35±0.25 


< 10 


812047 


0.2648 


3.26±0.57 


16.31±2.30 


< 10 


818329 


0.2656 


3.16±0.56 


12.90±1.83 


< 10 


825619 


0.3100 


2.30±0.49 


1.72±0.38 


< 10 


817135 


0.3220 


1.79±0.39 


4.48±0.71 


< 10 


822887 


0.2167 


4.56±0.92 


5.66±0.89 


< 10 


840778 


0.2614 


2.79±0.70 


1.84±0.55 


< 10 


848386 


0.2861 


2.29±0.56 


0.90±0.22 


< 10 


831576 


0.2512 


1.96±0.42 


0.94±0.21 


< 10 


841054 


0.2726 


1.91±0.48 


1.38±0.44 


< 10 


840771 


0.3302 


1.65±0.36 


0.83±0.17 


< 10 


836248 


0.2602 


2.66±0.57 


1.75±0.32 


< 10 


837327 


0.2190 


5.33±1.06 


4.03±0.66 


< 10 


826959 


0.2607 


1.73±0.39 


2.46±1.00 


< 10 


844388 


0.2158 


1.99±0.41 


4.02±0.61 


< 10 


814007 


0.5866 


1.75±0.75 


1.34±0.28 


10-11 


801253 


0.4803 


2.88±0.55 


2.28±0.36 


10-11 


812879 


0.2506 


7.73±1.35 


26.64±4.05 


10-11 


814176 


0.3753 


1.96±0.42 


1.38±0.28 


10-11 


813806 


0.6896 


1.65±0.67 


2.79±1.73 


10-11 


820021 


0.6753 


1.81±0.50 


0.52±0.27 


10-11 


826685 


0.6018 


2.42±0.56 


1.75±0.33 


10-11 


825318 


0.4789 


0.98±0.21 


1.28±0.28 


10-11 


817189 


0.4609 


3.10±0.69 


36.97±5.32 


10-11 


823751 


0.4795 


2.73±0.58 


2.54±0.43 


10-11 


827923 


0.4332 


1.37±0.30 


1.38±0.25 


10-11 


834174 


0.5027 


1.69±0.47 


2.33±0.51 


10-11 


833601 


0.3342 


2.49±0.55 


1.54±0.28 


10-11 


824223 


0.6776 


1.83±0.63 


1.58±0.69 


10-11 


832126 


0.3711 


3.18±0.62 


5.32±0.81 


10-11 


831349 


0.3115 


4.43±0.89 


3.56±0.56 


10-11 


835952 


0.3616 


2.31±0.52 


1.16±0.25 


10-11 


844837 


0.5021 


1.92±0.41 


1.97±0.32 


10-11 


843329 


0.5064 


2.56±0.46 


13.91±2.02 


10-11 


804831 


0.3466 


2.69±0.56 


1.32±0.24 


10-11 


834565 


0.5339 


1.82±0.54 


2.23±0.36 


10-11 


826188 


0.3453 


3.93±0.75 


3.71±0.65 


10-11 


816998* 


0.4248 


4.86±1.03 


8.62±1.36 


10-11 


838297 


0.5037 


1.22±0.36 


0.91±0.21 


10-11 


835862 


0.4027 


3.32±0.69 


1.78±0.34 


10-11 


844486 


0.4694 


1.93±0.40 


0.83±0.21 


10-11 


851898 


0.6655 


2.02±0.51 


4.57±0.83 


10-11 


826453 


0.3547 


2.76±0.56 


1.12±0.28 


10-11 


803488 


0.5867 


4.69±1.79 


39.11±9.32 


> 11 


820949 


0.5712 


3.00±0.60 


2.23±0.45 


> 11 


812432 


0.6611 


7.59±2.38 


39.35±11.11 


> 11 


817871 


0.6738 


2.42±0.53 


14.74±3.06 


> 11 


841690 


0.4821 


11.01±1.92 


58.98±8.40 


> 11 


831770 


0.6877 


3.35±0.61 


3.44±0.63 


> 11 


830349 


0.6693 


1.94±0.48 


1.59±0.42 


> 11 


838689* 


0.6884 


10.84±6.62 


24.96±5.88 


> 11 


836868 


0.6793 


2.42±0.61 


23.14±3.16 


> 11 


851888 


0.6667 


2.27±0.88 


3.07±0.97 


> 11 


820787 


0.6753 


3.25±0.58 


26.00±3.72 


> 11 
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TABLE 3 

The Ha A6563 and [Nil] A6584 line fluxes for the galaxies in our 

SAMPLE FOR WHICH BOTH LINE EW ARE GREATER THAN 5 A. ALL FLUXES ARE 
APERTURE AND EXTINCTION-CORRECTED. Ha A6563 FLUXES ALSO INCLUDE 
CORRECTIONS FOR STELLAR ABSORPTION. THE FLUX UNITS ARE 

10~ 15 crgs^crrL 2 . zCOSMOS BLAGN have been excluded. The SOURCE 

FLAGGED WITH AN ASTERISK CORRESPONDS TO AN X-RAY AGN WITH 
L x > 10 42 ' 5 ergs" 1 . 



zCOSMOS id zspec HaA6563 [Nil] A6584 log 10 uLZ' 
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Fig. 1 — The zCOSMOS-bright 10k and the SCOSMOS/MIPS shallow sample coverage fields. 
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Fig. 2. — Distributions of 24/xmfluxes of all the SCOSMOS-shallow 24 fim sources present in the zCOSMOS-bright lOk-sample field 
(empty histogram); those that have a counterpart in the zCOSMOS-bright parent catalogue (thin-line, shaded histogram); and those that 
have a good-quality redshift in the zCOSMOS-bright 10k catalogue (thick-line, shaded histogram). 
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Fig. 3. — Rest-frame 24 fim luminosities of the 609 galaxies in our zCOSMOS-lOk sample versus the spectroscopic redshifts. Red circles 
correspond to secure one-to-one associations and blue triangles indicate the cases in which there is more than one optical counterpart to 
the 24 lira source within 2 arcsec radius (with the zCOSMOS source being the closest to the 24 fim centroid). Crosses and plus-like symbols 
indicate spectroscopically classified BLAGN and IRAC power-law SED galaxies, respectively. 
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Fig. 4. — Composite average zCOSMOS spectra of 24 fim galaxies with vL^ A ' m < 10 10 Lq in different redshift bins. From top to bottom, 
the numbers of stacked galaxies are 74, 64 and 38, respectively. 
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FlG. 5. — Composite average zCOSMOS spectra of 24 ^im galaxies with 10 10 < uL^ ^ m < 10 11 Lq in different redshift bins. From top to 
bottom, the numbers of stacked galaxies are 70, 57 and 42, respectively. 
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FlG. 6. — Composite average zCOSMOS spectra of 24 fim galaxies with zaL^ 4 /J,m > 10 11 L0 in different rcdshift bins. From top to bottom, 
the numbers of stacked galaxies are 42, 86 and 89, respectively. 
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Fig. 7. — Rest-frame EW measurements of different lines in the zCOSMOS spectra of 24 fim galaxies as a function of redshift. Symbols 
of different size and colour correspond to galaxies with different mid-IR luminosities: !/L^ 4Mm < 10 10 L0 (small blue circles), 10 10 < 
vL v ' lm < 10 11 Lq (medium-size green circles) and vL v ' Jm > 10 11 Lq (large red circles). 
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Fig. 9. — Balmer decrement between the HaA6563and H/3 A4861 lines for our 24 fim galaxies at 0.2 < z < 0.3, versus the logarithm of 
their (far-IR+UV) over UV luminosity ratio. We consider the latter quantity as proportional to the extinction in the UV. Only galaxies 
with EW> 5 A for both the Ha A6563 and the H/3 A4861 lines are shown. The colour and size code for the circles is the same as in previous 
figures. Thick lines of different styles and colours indicate different reddenning laws, as detailed in the inset within the plot. The horizontal 
dashed line marks the value of the intrinsic Balmer decrement in a case B recombination with temperature T = 10, 000K. We note that 
some galaxies have L(Ha)/L(H/3) ratios that are equal to or lower than the case B recombination, but at the same time have significantly 
large (£70 /jm + £200o)/^2000 ratios. This is probably due to aperture effects, i.e. the f -arc sec wide V1MOS slit maps the unobscured 
central region of the galaxy, while the 1R flux is produced by surrounding dust (see Figure [TOt . 
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Fig. 10. — HST/ACS /-band postage stamps of six galaxies characterised by optical spectra with an Ha A6563 -to-H/3 A4861 ratio equal 
to or lower than the intrinsic absorption in the case B recombination. The size of each stamp is 7.5 X 7.5 arcsec. The position of the 
VIMOS slit is shown on each image. 




Fig. 11. — V-band extinctions derived for the 24 /an galaxies versus redshift, using the Calzetti et al. reddening law. In each case, we 
obtained the extinction by equating the SFR derived from the HaA6563or H/3 A4861 lines with the total (IR+UV) SFR. The colour and 
size code for the circles is the same as in previous figures. Triangles indicate galaxies from the CFRS (Lilly et al. 1995), for which SFR and 
extinctions have been computed by Maier et al. (2005). For the three CFRS galaxies that have been detected at 15 (im with ISO (Flores 
et al. 1999; filled triangles), we recomputed the SFR and extinctions taking into account their IR flux. 
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Fig. 12.— Left panel: [OIII] A5007 /H/3 A4861 ratios for our galaxies in the cases that both H/3A4861and [OIII] A5007 have EW > 5 A, 
versus redshift. The colour and size code for circles is the same as in previous figures. Filled circles with a plus sign within correspond to 
(narrow-line) X-ray AGN with L x > 10 42 5 erg s" 1 . Measurements for those cases with H/3A4861EW > 5 A but [OIII] A5007EW < 5 A 
are indicated with small empty circles. Asterisks correspond to our measurements on average stacked spectra. The colour and size code 
for the asterisks is equivalent to that for circles. Brown dots correspond to measurements performed on all other zCOSMOS galaxies with 
H/3 A4861 EW > 5 A at similar redshifts (Lamareille et al., in preparation). Small green and large red crosses correspond to measurements 
listed in Table 2 of Rupke et al. (2008) for LIRGs and ULIRGs, respectively. Right panel: [OIII] A5007 /H/3 A4861 ratios for sub-millimetre 
galaxies at z > 1.5 obtained by Takata et al. (2006). 
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Fig. 13.— [Nil] A6584 /Ha A6563 line ratios versus redshift for our galaxies in the cases that both Ha A6563 and [Nil] A6584 have EW > 5 A. 
The colour and size code for circles is the same as in previous figures. Small empty circles indicate measurements for sources with 
Ha A6563 EW > 5 A but [Nil] A6584 EW < 5 A. 
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Fig. 14.— The BPT diagram for our 24 fim sources at 0.2 < z < 0.3 with EW > 5 A for the four H/3A4861 , [OIII] A5007 , H«A6563and 
[Nil] A6584 lines. Small empty circles indicate sources for which only HaA6563and H/3 A4861 have EW > 5 A. Asterisks indicate measure- 
ments on composite spectra at median redshifts z = 0.25 and 0.35. The colour and size code for circles and asterisks is the same as in 
previous figures. Brown dots correspond to measurements performed on all other zCOSMOS galaxies with HoA6563and H/3A4861EW 
> 5 A at similar redshifts (Lamareille et al., in preparation). The dashed line indicates the starburst/AGN division line proposed by Kewley 
et al. (2001), while the dot-dashed line corresponds to the limit for composite systems (Kauffmann et al. 2003a). 
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FlG. 15. — The same composite spectra shown in the bottom panel of figure [4] and in the top panel of figure [5] but in this case extended 
beyond observed wavelength A = 8500 A, i.e. the region significantly affected by fringing. The composite has been made by stacking the 
spectra of 38 and 70 galaxies, respectively, with i/L^ 4Alm < 10 10 Lq and 10 10 < z^L^ 4AIm < 10 11 Lq , and with median redshifts z = 0.34 
and 0.35. 
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Fig. 16. — Oxygen abundances derived for 24 /xm/zCOSMOS galaxies at different redshifts. All the abundances have been calculated 
with the same technique, assuming the dust extinction to be known and fixed (see text). Cases of non-degenerate abundances are indicated 
by filled circles with errors bars. Upward-pointing arrow indicate lower limits. Cases of degenerate metallicities are represented by empty 
circles (upper branch) and crosses (for the corresponding lower- branch values). The colour and size code for circles is the same as in 
previous figures. The circle with a plus sign within corresponds to a (narrow-line) X-ray AGN with Lx > 10 42 ' 5 ergs -1 . The dashed line 
indicates the value of the solar metallicity determined by Asplund et al. (2004). 
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Fig. 17.— [OIII] A5007 /H/3 A4861 versus [Oil] A3727 /H/3 A4861 ratios for the galaxies with derived metallicities at 0.5 < z < 0.7 shown in 
figures ri6l and [T8l Filled circles indicate cases of non-degenerate metallicities, while empty circles correspond to cases with degenerate values. 
The solid and dashed lines are, respectively, the median and confidence limits of the empirical separation between star-forming galaxies 
and AGN derived by Lamareille et al. (2004). The (narrow-line) X-ray AGN with L x > 10 42 5 ergs -1 has logio([OII/H/3]) = 0.81 ± 0.20 
and logio([OIII/H/3]) = 0.36 ± 0.20 (red large circle behind the green small one). 
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Fig. 18. — Oxygen abundances versus already-assembled stellar masses for our 24 fim -selected galaxies. Top panel: 0.2 < z < 0.3 
galaxies. Bottom panel: 0.5 < z < 0.7 galaxies. The symbol code is the same as in previous figures. Triangles in the bottom panel indicate 
galaxies from the CFRS (Lilly et al. 1995) with known metallicities (Maier et al. 2005). Those triangles filled in yellow correspond to 
CFRS galaxies with ISO 15 /an detections (Flores et al. 1999). The solid line shows the best-fit mass-metallicity relation obtained by 
Tremonti et al. (2004) on SDSS data. The dashed lines correspond to its evolution at 0.5 < z < 0.7 and 0.7 < z < 1.0, as derived from all 
zCOSMOS-bright 10k galaxies (Contini et al. 2008). 
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Fig. 19. — Left panel: redshift distributions of all BLAGN in the zCOSMOS-bright 10k catalogue (empty histogram) and those detected 
in the MlPS-shallow catalogue (filled histogram). Right panel: 24 /im luminosity distributions of all zCOSMOS-bright BLAGN in our 
sample (thick-line, empty histogram) and the narrow line IR sources associated with X-ray AGN with Lx > 10 42 ' 5 ergs -1 (thin-line, filled 
histogram). 
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Fig. 20. — Composite spectrum of 8 IRAC power-law SED 24 /imsources at 0.5 < z < 0.85 whose individual spectra are not classified as 
BLAGN. 
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Fig. 21. — The expected evolution of H<5 versus D n (4000) for galaxies with different SFH: exponentially-declining (with r = 0.01 and 
0.1 Gyr; solid and dotted lines, respectively) and a constant star formation (dashed line). The evolution has been measured using the 
Bruzual & Chariot (Bruzual 2007) synthetic templates. 
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Fig. 22. — D n (4000) values for 24 (im galaxies as a function of redshift. Brown dots indicate all the other zCOSMOS galaxies in the same 
redshift range (Lamareille et al., in preparation). 
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Fig. 23. — The average location of 24 fim galaxies in the H<5 versus D n (4000) diagram. The c olou r and size code for the asterisks is the 
same as in previous figures. Thick lines represent the same modelled evolution shown in figure 1211 The thin solid lines show the effect of 
secondary bursts of star formation that produce 5 or 10% of the already-assembled stellar mass. We show examples where the burst starts 
when the galaxy is 2 or 3 Gyrs old. The numbers indicated with a plus sign in front correspond to the time passed since the onset of the 
secondary burst. 




Fig. 24. — Example of a 24 /an galaxy candidate to be at the earliest stages of a burst of star formation. Left panel: spectrum over the 
total observed wavelength range. Right panel: detailed view of the region of high-order Balmer lines. 



